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These introductory lecture notes describe recent results on novel Mott insulating phases which 
are "descendents" of superconductors - obtained by "quantum disordering" . After a brief overview 
of quantum magnetism, attention is focussed on the spin-liquid phase of the two-leg Hubbard ladder 
and the nodal liquid - a descendent of the d x 2_ y i superconductor. These notes, which will appear 
in a future Les Houches publication, are self-contained and an effort has been made to keep them 
accessible. 



I. INTRODUCTION 

At the foundation of the quantum theory of metals is 
the theory of the non-interacting electron gas, in which 
the electrons move through the material interacting only 
with the periodic potential of the ions, and not with 
one another. Surprisingly, the properties of most met- 
als are quite well described by simply ignoring the strong 
Coulomb repulsion between electrons, essentially because 
Pauli exlusion severely limits the phase space for electron 
collisions. EJ But in some cases electron interactions can 
have dramatic effects leading to a complete breakdown of 
the metallic state, even when the conduction band is onbz 
partially occupied. In the simplest such Mott insulatorp 
there is only one electron per crystalline unit cell, and so 
a half-filled metallic conduction band would be expected. 

With the discovery of the cuprate superconductors in 
1 986,0 there has been a resurgence of interest in Mott 
insulators. There are two broad classes of Mott insu- 
lators, distinguished by the presence or absence of mag- 
netic order QU More commonly spin rotational invariance 
is spontaneously broken, and long-range-jmagnetic order, 
typically antiferromagnetic, is realized.EJ There are then 
low energy spin excitations, .-the spin one magnons. Al- 
ternatively, in a spin-liquidB Mott insulator there are 
no broken symmetries. Typically, the magnetic order is 
short-ranged and there is a gap to all spin excitations : 
a spin-gap. 

In the cuprates the Mott insulator is antiferromagnet- 
ically ordered, Lm but upon doping with holes the anti- 
ferromagnetism is rapidly destroyed, and above a cer- 
tain level superconductivity occurs with d x i_ y i pairing 
symmetry. But at intermediate doping levels between 
the magnetic and d-wave superconducting phases, there 
are experimental signs of a spin gap opening below a 
crossover temperature T*(x) (see Figure 1). The ulti- 
mate nature of the underlying quantum ground state in 
this portion of the phase diagram - commonly called the 
pseudo-gap regime - is an intriguing puzzle. More gener- 
ally, the apparent connection between a spin-gap and su- 
perconductivity has been a source of motivation to search 
for Mott insulators of the spin-liquid variety. 



Generally, spin liquids are more common in low di- 
mensions where quantum fluctuations can suppresa_mag- 
netism. Quasi-one-dimensional ladder materialsoEj are 
promising in this regard and have received extensive at- 
tention, particularly the two-leg ladder.El The Mott in- 
sulating spin-liquid phase of the two-leg ladder can be 
understood by mapping to an appropriate spin-model - 
the Heisenberg antiferromagnet. Spin- liquid behavior re- 
sults from the formation-oLsinglct bond formation across 
the rungs of the ladder .Effij 

Almost without exception, theoretical studies of spin- 
liquids start by mapping to an appropriate spin-model, 
and the charge degrees of freedom are thereafter ignored. 
This represents an enormous simplification, since spin 
models are so much easier to analyze that the underlying 
interacting electron model. This approach to quantum 
magnetism has yielded tremendous progress in the past 
decadeB But is the simplification to a spin-model always 
legitimate? A central goal of these lectures is to ana- 
lyze a novelptwei-dimensional spin-liquid phase - called a 
nodal liquidr ti 15 i - which cannot be described in terms of 
a spin model. Although the nodal liquid is a Mott insu- 
lator with a charge gap and has no broken symmetries, 
it possesses gapless Fermionic degrees of freedom which 
carry spin. 
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Fig. 1: Schematic phase diagram of a high-temperature 
superconductor as a function of doping x and tempera- 
ture T. 
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The standard route to the spin-liquid invokes quantum 
fluctuations-,to suppress the magnetic order of a quantum 
spin-modeLu The proximity of antiferromagnetism to d- 
wave superconductivity in the cuprates suggests an alter- 
nate route. Indeed, as we shall see, the nodal liquid phase 
results when a d — wave superconductor is "quantum dis- 
ordered" . The gapless Fermionic excitations in the nodal 
liquid are descendents of the low energy quasiparticles of 
the d-wave superconductor. 

The spin-liquid phase of the two-leg ladder gives us 
a simpler example of a quantum disordered supercon- 
ductor. To see this, we will revisit the two-leg lad- 
der, employing a model of interacting electrons^ rather 
than truncating to a spin-model. Retaining the charge 
degrees of freedom will enable us to show that the 
Mott-insulating phase of the two-leg ladder actually 
exhibits pairing, with an approximate d-wave symme- 
try. Moreover, upon doping, the two-leg ladder exhibits 
quasi-long-range superconducting (d-wave) pairing cor- 
relations. This behavior is reminiscent of that seen in 
the underdoped cuprate superconductors. 

These notes are organized as follows. In Section |H] 
a simple tight binding model of interacting electrons is 
introduced and it's symmetry properties are discussed. 
Section [II specializes to the Mott insulating state at 



half-filling, focussing on the magnetic properties employ- 
ing the Heisenberg antiferromagnet spin-model. In Sec- 
tion |^ the method of Bosonization is briefly reviewed 
for the case of a one-dimensionless spinless electron gas. 
Section [v] is devoted to an analysis of the Mott insulating 
state of the two-leg Hubbard ladder, employing a weak 
coupling perturbative renormalization group approach. 
The remaining sections focus on the two-dimensional d- 
wave superconductor, and the nodal liquid phase which 
descends from it upon quantum disordering. Specifically, 
Section VI briefly reviews BCS theory for a d-wave su- 



perconductor focussing on the gapless quasiparticles. An 
effective field theory for this state, includ ing quantum 



phase fluctuations, is obtained in Section VII. Section 



VIII implements a duality transformation of this effec- 



tive field theory, which enables a convenient description 
of the nodal liquid phase in Section IX. 



band are built up from a single atomic orbital on each 
of the ions in the solid. The conduction electrons are 
presumed to be moving through the solid, tunnelling be- 
tween ions. We denote the creation and annihilation op- 
erators for an electron with spin a =|,| on the ion at 
position x by c\{x) and c a (x). These operators satisfy 
the canonical Fermionic anti-commutation relations, 



Ca(x),ci(x')}_ = SapSx 



(2.1) 



If the orbitals in question form a simple Bravais lattice 
with, say, cubic symmetry, then the appropriate tight 
binding Hamiltonian is, 

H = -tJ2 [4(s)cjs0 + kc.]-/*5>(a0, (2.2) 

(xx') x 

where the first summation is over near neighbor sites. 
Here t is the tunnelling rate between neighboring ions 
and for simplicity we have ignored further neighbor tun- 
nelling. The electron density n(x) = c] x (x)c a {x) can be 
adjusted by tuning the chemical potential, /x. 

In the Cuprate superconductors Copper and Oxygen 
atoms form two dimensional sheetsJH with the Copper 
atoms sitting at the sites of a square lattice and the Oxy- 
gen atoms sitting on the bonds, as depicted schematically 
in Figure 2. In the simplest one-band models the sites of 
the tight binding model are taken as the Copper atoms, 
and c' (x) removes an electron (adds a hole) from a Cop- 
per 3d orbital. In most of the materials the 3d shell is 
almost filled with roughly one hole per Copper atom, so 
that the tight binding model is close to half-filling with 
(n(x)) « 1. 

The tight binding Hamiltonian is invariant under 
translations by an arbitrary real space lattice vector, R, 



(x) 



R). 



(2.3) 



This discrete symmetry implies the conservation of crys- 
tal momentum, up to a reciprocal lattice vector, G, with 
exp{iG ■ R) = 1. Being quadratic, the Hamiltonian can 
be diagonalized by transforming to (crystal) momentum 
space by defining, 



II. MODELS AND METALS 



c a (x) 
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(2.4) 



A. Noninteracting electrons 

In metals the highest lying band of Bloch states is only 
partially occupied, and there are low energy electronic 
excitations which consist of exciting electrons from just 
below the Fermi energy into unoccupied states. These 
excitations can be thermally excited and contribute to 
thermodynamic properties such_as the specific heat, as 
well as to electrical conduction.!!] Tight binding models 
give a particularly simple description of the conduction 
band. In the simplest case the states in the conduction 



Here V denotes the "volume" of the system, equal to the 
total number of sites N with the lattice spacing set to 
unity, and the sum is over crystal momentum within the 
first Brillouin zone compatible with periodic boundary 
conditions. The momentum space creation and anihilla- 
tion operators also satisfy canonical Fermion anticommu- 
tation relations: 



(2.5) 



In momentum space the Hamiltonian takes the standard 
diagonal form, 
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efcC fca Cfco 



(2.6) 



invariant under the discrete translation symmetry: 
Cfect — > e tk R Cka- For a 2d square lattice with near- 
neighbor hopping, the energy is simply 



efc = — 2t[cosk x + cos k y ] — /i. 



(2.7) 



The ground state consists of filling those states in mo- 
mentum space with negative, leaving the positive en- 
ergy states unoccupied. The Fermi surface, separating 
the occupied from empty states, is determined by the 
condition ej- = 0. For the 2d squa re lattice at half-filling 
with energy dispersion Eqn. 2.7 (at \i = 0), the Fermi 
surface is a diamond, as shown in Figure 2. 




Fig. 2: Schematic illustration of a single Copper- 
Oxygen plane, consisting of a square lattice of Copper 
atoms (solid points) and Oxygen atoms (open circles). 
Two-dimensional Brillouin zone for the 2d square lat- 
tice tight binding model with near neighbor hopping is 
shown at right. At half- filling all states in the Fermi sea 
(shaded) are occupied. 

Particle/hole excitations above the ground state con- 
sist of removing an electron from within the full Fermi 
sea, and placing it in an unoccupied positive energy state. 
In most metals the width of the conduction band (propor- 
tional to t) is of order an electron volt (roughly 10 4 K) so 
that even at room temperature only "low energy" par- 
ticle/hole states confined within close proximity to the 
Fermi surface are thermally excited. In addition to be- 
ing thermally active, these low energy particlc/holc ex- 
citations can be excited by an electric field, and lead to 
metallic electrical conduction. 

In the band theory of solids, insulators occur whenever 
the highest lying energy band is fully occupied. Excited 
states then involve promoting electrons into the next 
available band which typically requires a very large en- 
ergy (electron volts). Not surprisingly, such band insula- 
tors are very poor conductors of electricity. By constrast, 
in Mott insulators the highest band is only partially oc- 
cupied, yet conduction is blocked by strong electron in- 
teractions. 

Before addressing the complications of electron interac- 
tions, it is instructive to briefly consider the symmetries 
of the above Hamiltonian, and the associated conserved 
quantities. There are only two continous symmetries, 



associated with conservation of charge and spin. The 
Hamiltonian is invariant under the global U(l) charge 
symmetry, 



c a {x) 



3 c Q (32) , 



(2.8) 



for arbitrary (constant) angle 9q. Conservation of spin is 
due to the global SU{2) symmetry, c a {x) — > U a pc^(x), 
with U = cxp(i0-er) and Pauli matrices cr a j3. The Hamil- 
tonian is invariant under this transformation, Hq — > Ho, 
for arbitrary spin rotations 6. Here and below we ig- 
nore spin-orbit effects which (usually weakly) break the 
continuous spin rotational symmetry. 

There are also a number of discrete symmetries. The 
Hamiltonian is real, Hq = Hq, a signature of time re- 
versal invariance (for models with spin-independent in- 
teractions). For a square lattice the Hamiltonian is also 
invariant under reflection (or parity) symmetry, c a {x) — > 
c a (—x). This implies that gfc = On the square 

lattice, a discrete particle/hole transformation is imple- 
mented by 



t (x) 



p/h 



<&(*), 



(2.9) 



with re — (jr, 7r). At half-filling when /i = 0, H is in- 
variant under this symmetry, but with further neighbor 
hopping terms the kinetic energy will generally not be 
particle/hole symmetric. In momentum space the parti- 
cle/hole transformation is implemented via Cka - > c lr-ka 
and invariance of the kinetic energy implies that e?c = 



B. Interaction Effects 

Spin-independent density interactions can be included 
by adding an additional term to the Hamiltonian: 



H\ = - £^ v(x — x')n{x)n(x'). 



(2.10) 



x,x' 



For Coulomb interactions v(x) ~ e 2 /|x| is long-ranged. 
For simplicity the long-rangecLiiiteractions are often ig- 
nored. In the Hubbard modeloa only the on-site repul- 
sive interaction is retained, 



H u = u^2rif(x)ni(x), 



(2.11) 



with 



= A 



(y a c a . This can be re-cast into a manifestly 



spin-rotationally invariant form: 



(2.12) 



Despite the deceptive simplicity of these effective models, 
they are exceedingly difficult to analyze. Even the Hub- 
bard Hamiltonian, H — Hq + H u , which is parameterized 
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by just two energy scales, t and u, is largely intractable,cl 
except in one-dimension. Since the typical interaction 
scale u is comparable to the kinetic energy t there is no 
small parameter. Moreover, one is typically interested 
in phenomena occuring on temperature scales which are 
much smaller than both u and t. 

In most metals, the low energy properties are quite 
well described by simply ignoring the (strong!) interac- 
tions. This surprising fact can be understood (to some 
degree) from Landau's Fermi-liquid thcotyjj and more re- 
cent renormalization group arguments.tll The key point 
is that the phase space available for collisions between 
excited particles and holes vanishes with their energy. 
In metals the phase space is evidently so restrictive that 
the surviving interactions do not change the qualitative 
behavior of the low energy particle/hole excitations. In- 
deed, the quasiparticle excitations within Landau's Fermi 
liquid theory have the same quantum numbers as the 
electron (charge e spin 1/2 and momentum), but move 
with a "renormalized" velocity. But some materials such 
as the Cuprates are not metallic, even when band struc- 
ture considerations would suggest a partially occupied 
conduction band. In these Mott insulators one must in- 
voke electron interactions. 



III. MOTT INSULATORS AND QUANTUM 
MAGNETISM 

The Hubbard model at half-filling is perhaps the sim- 
plest example of a Mott insulator. To see this, consider 
the behavior as the ratio u/t is varied. As discussed 
above, for u/t — the model is diagonalized in momen- 
tum space, and exhibits a Fermi surface. But at half- 
filling the model is also soluble when u/t = oo. Since 
the onsite Hubbard energy takes the form, u(n — l) 2 /2, 
in this limit the ground state consists simply of one elec- 
tron on each site. The electrons are frozen and immobile, 
since doubly occupied and unoccupied sites cost an en- 
ergy proportional to u. The state is clearly insulating - 
a Mott insulator. 

In this large u limit it is very costly in energy to add 
an electron, and the state exhibits a charge gap of order 
u. But there are many low energy spin excitations, which 
consist of flipping the spin of an electron on a given site. 
For infinite u this spin-one excitation costs no energy 
at all, and indeed the ground state is highly degenerate 
since the spins of each of the N localized electrons can 
be either up or down. For large but finite u/t one still 
expects a charge gap, but the huge spin degeneracy will 
be lifted. 

The fate of the spin degrees of freedom in the Mott 
insulator is enormously interesting. Broadly speaking, 
Mott insulators come in two classes, distinguished by the 
presence or absence of spontaneously broken symmetries. 
Often the spin rotational invariance is spontaneously bro- 
ken and the ground state is magnetic, but SU(2) invari- 
ant spin structures which break translational symmetries 



are also possible. In the second class, usually referred to 
as spin liquid states there are no broken symmetries. 

A. Spin Models and Quantum Magnetism 

Traditionally, spin physics in the Mott insulating states 
have been analyzed by studying simple spin models. 
These focus on the electron spin operators: 

s ( x ) = T^cii^vapcpix), (3.1) 

where er is a vector of Pauli matrices. These spin oper- 
ators satisfy standard angular momentum commutation 
relations: 

[S tl (x),S„(x']=i5 xx ,e^ x S x . (3.2) 
They also satisfy, 

S 2 {x) = |n(x)[2-n(x)]. (3.3) 

Within the restricted sector of the full Hilbert space with 
exactly one electron per site, these operators are bone fide 
spin 1/2 operators satisfying S 2 = s(s+ 1) with s = 1/2. 
Their matrix elements in the restricted Hilbert space are 
identical to the Pauli matrices: <r/2. 

The simplest spin model consists of a (square) lattice of 
spin 1 /2 operators coupled via a near neighbor exchange 
interaction, J: 

H = J^2 S{x)-S(x'). (3.4) 

{xx') 

This spin model can be obtained from the half-filled Hub- 
bard model,[j by working perturbatively in small t/u. For 
t/u = the spins are decoupled, but an antifcrromagnetic 
exchange interaction J = 4t 2 /u is generated at second 
order in t. Specifically, the matrix elements of the spin 
Hamiltonian in the restricted Hilbert space arc obtained 
by using second order perturbation theory in t. The in- 
termediate virtual states are doubly occupied, giving an 
energy denominator u. 

Mapping the Hubbard model to a spin model repre- 
sents an enormous simplification. The complications due 
to the Fermi statistics of the underlying electrons have 
been subsumed into an exchange interaction. The spin 
operators are essentially bosonic, commuting at different 
sites. It should be emphasized that at higher order in t/u 
multi-spin exchange interactions will be generated, also 
between further separated spins. If t/u is of order one, 
then it is by no means obvious that it is legitimate to 
truncate to a spin model at all. 

A central focus of quantum magnetism during the past 
decade has been exploring the possible ground states and 
low energy excitations of such spin models .□'Q The above 
s = 1/2 square lattice Heisenberg antiferromagnet is, of 
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course, only one member of a huge class of such mod- 
els. These models can be generalized to larger spin s, 
to different lattices and/or dimensionalities, to include 
competing or frustrating interactions, to include multi- 
spin interactions, to "spins" in different groups such as 
SU(N), etc.. Not surprisingly, there is an almost equally 
rich set of possible ground states. 

The main focus of these notes is the 2d "nodal liq- 
uid" , a spin-liquid phase obtained by quantum disorder- 
ing a d-wave superconductor. As we shall see in Sec. IX, 
in the nodal liquid the spin excitations are carried by 
Fermionic degrees of freedom and cannot be described by 
(Bosonic) spin operators. In truncating to the restricted 
Hilbert space with one electron per site, one has effec- 
tively "thrown out the baby with the bath water" . The 
nodal liquid phase requires retaining the charge degrees 
of freedom. 

But spin models are much simpler than interacting 
electron models, relevant to many if not most Mott in- 
sulators (as well as other localized spin systems) and 
extremely rich and interesting in their own right. So 
I would like to briefly summarize some of the possi- 
ble ground statesy-|focussing on spin 1/2 models on bi- 
partite lattices. BO Consider first those ground states 
with spontaneously broken symmetries. Most common 
is the breaking of spin-rotational invariance. If the spin 
operators are treated as classical fixed length vectors, 
which is valid in the large spin limit (s — > oo), the ground 
state of the near neighbor square lattice antiferromagnet 
is the Neel state (up on one sublattice, down on the other) 
which breaks the SU(2) symmetry. For finite s the Neel 
state is not the exact ground state, but the ground state 
is still antiferromagnetically ordered, even for s = 1/2. 
Quantum fluctuations play a role in reducing the sub- 
lattice magnetization, but (for the 2d square lattice) do 
not drive it to zero. The low energy excitations are gap- 
less spin- waves (ie. magnons), as expected when a con- 
tinuous symmetry is spontaneously broken. 

For some spin models the ground state is spin ro- 
tationally invariant but spontaneously breaks (discrete) 
translational symmetry. The n classic example is the 
Majumdar-Ghosh Hamiltonian,u 

Hmg = J^ S ( X ) ■ S ( x + + \ S{ - X ) ■ S ( x + 2 )]> 

X ~ 

(3.5) 

which describes a one dimensional s — 1/2 Heisenberg 
antiferromagnetic spin chain with a second neighbor ex- 
change interaction. The exact ground state of this model 
is a two-fold degenerate "spin-Peierls" state: 

|G)± = J]!! Tto)| 42.il) - | i2x)| Ta*±i>]. (3.6) 



This state consists of a product of "singlet bonds" formed 
from neighboring pairs of spins, and breaks invariance 



under translations by one lattice spacing. Since the sin- 
glet bonds are rotationally invariant, the SU(2) symme- 
try remains unbroken. The second neighbor interaction 
has effectively suppressed the tendency towards antifer- 
romagnetic order. 



B. Spin Liquids 

Spin liquid ground states in which no symmetries are 
broken generally occur more readily in low dimensions 
where quantum fluctuations are more effective at de- 
stroying magnetic order. The one-dimensional s = 1/2 
chain with near neighbor antiferromagnetic exchange ex- 
hibits power law magnetic correlations at the antiferro- 
magnetic wave vector 7t.e3 Although "almost" magneti- 
cally ordered the SU (2) symmetry is not broken in the 
ground state, which thus technically qualifies as a spin 
liquid. More dramatic is the behavior of the s = 1/2 an- 
tiferromagnetic two-leg ladder, shown in Figure 3. This 
model exhibits a featureless spin-rotationally invariant 
ground state with exponentially decaying spin correla- 
tion fuctioips.. and a non-zero energy gap for all spin 
excitations.E£lE3 The physics can be best understood in 
the limit in which the exchange interaction across the 
rungs of the ladder greatly exceeds the intra-leg ex- 
change: Jj_ >> J. When J = the ground state consists 
of singlet bonds formed across the rungs of the ladder, 
with triplet excitations separated by an energy gap of 
order J±. Perturbing in small J will cause these singlet 
bonds to "resonate" , but one expects the spin gap to sur- 
vive at least for J << J±. It turns out that the ground 
state evolves adiabatically and smoothly with increasing 
J, and in fact the spin-liquid survives for arbitrarily large 
J±/J. 



Ji. 



Fig. 3: Heisenberg spin model on a two-leg ladder. 
Spin 1/2 operators sit on the sites of the ladder, interact- 
ing via an antiferromagnetic exchange J along the ladder 
and J_l across the rungs. 

There has been an enormous amount of theoretical 
effort expended searching for two-dimensional spin 1/2 
models which exhibit spin-liquid ground states analogous 
to the two-leg ladder - but with little success. The orig- 
inal motivation soon after the discovery of superconduc^ 
tivity in the Cuprates was based on Anderson's ideascil 
that a Mott insulating spin-liquid exhibits "pre-formed" 
Cooper pairing. Doping the Mott insulator would give 
the Cooper pairs room to move and to condense into a 
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superconducting state, presumed to have s-wave pairing 
symmetry. But it soon became clear that the undoped 
Mott insulator in the Cuprates is not a spin-liquid, but 
actually antifcrromagnetically ordered. Moreover, recent 
experiments have established that the pairing symmetry 
in th e sup erconducting phase is d-wave rather than s- 
wave.liaE] 

However, recent theoretical workllS (see Section [v| be- 
low) has established that the pairing in the spin-liquid 
phase of the two-leg ladder actually has (approximate) 
d-wave symmetry. Moreover, doping this Mott,.iiisula- 
tor does indeed give the pairs room to move , r HP 2 ! and 
they form a one-dimensional d-wave "superconductor" 
(with quasi- longj ranged pairing correlations). The nodal 
liquid phaseoiij discussed extensively below is a two- 
dimensional analog of this spin-liquid phase. Indeed, 
we shall explicitly construct the nodal liquid by quan- 
tum disordering a two-dimensional d-wave superconduc- 
tor. As we shall see, the resulting 2d nodal liquid posesses 
gapless Fermionic excitations, which are descendents of 
the d-wave quasiparticles. These Fermions carry spin but 
no charge. The nodal liquid presumably cannot be the 
ground state of any (Bosonic) spin-model. To describe 
the nodal liquid one must employ the underlying inter- 
acting electron model which retains the charge degrees of 
freedom. 

Recent experiment has focussed attention on the un- 
derdoped regime of the Cuprate materialsjs occuring be- 
tween the antiferromagnetic and superconducting phases 
(see Fig. 1). In this pseudo-gap regime insulating behav- 
ior is seen at low temperatures, and there are indications 
for a spin gap - behavior reminiscent of a Mott insulating 
spin-liquid. We have suggestedliJ that this strange phase 
can perhaps be understood in terms of a doped nodal 
liquid. 

Before discussing further the 2d nodal liquid, it is in- 
structive to revisit the spin liquid phase of the two-leg 
ladder and analyze it directly with a model of inter- 
acting electrons. Specifically, we consider weak interac- 
tions (small u/t), a limit in which truncation to a spin 
model is not possible. This analysis is greatly aided by 
"Bosonization" - a powerful method which enables an 
interacting electron model in one dimension to be re- 
formulated in terms of collective Bosonic degrees of free- 
dom. See Ref . ^fj] as well as Fradkin's booka for useful 
reviews of Bosonization. First, in Sec. IV we briefly re- 



view Bosonization for the simplest case of a spinless one- 
dimensional electron gas, before turning to the two-leg 
ladder in Sec. 0. 



IV. BOSONIZATION PRIMER 



with hopping strength t. One can diagonalize this Hamil- 
tonia n by Fourier transforming to momentum space as in 
Eqn. 



2.4 



giving 



H 



k 



(4.2) 



with energy dispersion = —tcos(k) for momentum 
|fc| < 7r, as shown in Figure 4. In the ground state all of 
the negative energy states with momentum |fc| < k F are 
occupied. At half-filling the Fermi wavevector kp = ir/2. 
An effective low energy theory for these excitations can 
be obtained by focussing on momenta close to ±fcp- and 
defining continuum Fermi fields: 



i> R {q) = Ck 



F+qi 



^(?) = C- kF+q . 



(4.3) 



Here the subscripts R/L refer to the right/left Fermi 
points, and q is assumed to be smaller than a momen- 
tum cutoff, \q\ < A with A << kp. One can then 
linearize the dispersion about the Fermi points, writing 
( -±k F +q = ±vpq with vf the Fermi velocity. It is conve- 
nient to transform back to real space, defining fields 



M*) = -j= E eiqx M<i), 



(4.4) 



(with P = R,L) which vary slowly on the scale of the lat- 
tice spacing. This is equivalent to expanding the lattice 
electron operators in terms of continuum fields, 



c{x) ~ ip R (x)e ikFX + 4> L {x)e- tkFX . 



(4.5) 



After linearization, the effective low energy Hamiltonian 
takes the form, H = J dxTC, with Hamiltonian density, 



H = -v F [^ R id x il) R - iplid x i/j L 



(4.6) 



describing a one-dimensional relativistic Dirac particle. 
The associated Lagrangian density is simply 



C = ip R id t ip R + ij>lid x ip t 



H. 



(4.7) 




Consider the Hamiltonian for non-interacting spinless 
electrons hopping on a Id lattice, 



H 



-t^2^(x)c(x + 1) + h.c. 



(4.1) 



Fig. 4: Schematic energy dispersion for the one- 
dimensional electron gas. The negative energy states are 
occupied, with momentum |fc| < kp. The dispersion can 
be linearized around ±fcp, leading to a continuum Dirac 
Fermion theory. 
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5 ' We follow the heuristic development of 
which reveals the important physics, dispens- 



Consider a particle/hole excitation about the right 
Fermi point, where an electron is removed from a state 
with k < kp and placed into an unoccupied state with 
k + q > kp. For small momentum change q, the en- 
ergy of this excitation is u q = vpq. Together with 
the negative momentum excitations about the left Fermi 
point, this linear dispersion relation is identical to that 
for phonons in one-dimension. The method of Bosoniza- 
tion exploits this similarity by introducing a phonon dis- 
placement, fkld , 9, to decribe this linearly dispersing den- 
sity wa' 
Haldam 

ing with mathematical rigor. n To this end, consider a 
Jordan- Wigner transformationQ which replaces the elec- 
tron operator, c(x), by a (hard-core) boson operator, 

c(x) = 0{x)b{x) = exp[iTT n{x')]b(x). (4.8) 

x' <x 

where n(x) = c*{x)c(x) is the number operator. One 
can easily verify that the Bose operators commute at dif- 
ferent sites. Moreover, the lattice Hamiltonian Eqn. 4.1 
can be re-expressed in terms of these Bosons, and takes 
the identical form with c's replaced by b's. This trans- 
formation, exchanging Fermions for Bosons, is a special 
feature of one-dimension. The Boson operators can be 
(approximately) decomposed in terms of an amplitude 
and a phase, 



b(x) -»■ y/pe 



lip 



(4.9) 



We now imagine passing to the continuum limit, fo- 
cussing on scales long compared to the lattice spac- 
ing. In this limit we decompose the total density as, 
p(x) = po + p, where the mean density, po = kp/n, and p 
is an operator measuring fluctuations in the density. As 
usual, the density and phase are canonically conjugate 
quantum variables, taken to satisfy 



[cp(x),p(x')] = iS(x - x'). 



(4.10) 



Now we introduce a phonon-like displacement field, 9{x), 
via p(x) = d x 9(x)/ir. The full density takes the form: 
irp(x) = kp + d x 9. The above commutation relations are 
satisfied if one takes, 



[<p(x),6(x')] = -ine(x' - x). 



(4.11) 



Here O(x) denotes the heavyside step function, not to 
be confused with the displacement field 9. Notice that 
d x (p/iT is the momentum conjugate to 9. 

The effective (Bosonized) Hamiltonian density which 
describes the Id density wave takes the form: 



g- 1 {d x 9f 



(4.12) 



This Hamiltonian describes a wave propagating at ve- 
locity v, as can be readily verified upon using the com- 
mutation relations to obtain the equations of motion, 



d\ :9 = v d 2 9, and similarly for <p. Clearly one should 
equate v with the Fermi velocity, vp. The additional di- 
mensionless parameter, g, can be determined as follows. 
A small variation in density, p, will lead to a change in 
energy, E — p 2 /2k, where n = dp/ dp, is the compressibil- 
ity. Since d x 9 = irp, one deduces from TL that k = g/irv. 
But for a non- interacting electron gas, irvn = 1, so that 
9=1. In the presence of (short-ranged) interactions 
between the (spinless) electrons, one can argue that the 
above Hamiltonian density remains valid, but with renor- 
malized values of both g and v. This Hamilt»HK«i would 
then describe a (spinless) Luttinger liquid,E3'E3 rather 
than the free electron gas. 

The power of Bosonization relies on the ability to re- 
express the electron operator c(x) in terms of the Boson 
fields. Clearly c(x) must remove a unit charge (e) at x, 
and satisfy Fermion anticommutation relations. Consider 
first the Bose operator, b ~ exp(i^), which removes unit 
charge. To see this, note that one can write, 



(4.13) 



where P = d x <p/Tr is the momentum conjugate to 9. Since 
the momentum operator is the generator of translations 
(in 9), this creates a kink in 9 of height it centered at 
position x - which corresponds to a localized unit of 
charge since the density p = d x 9/w. To construct the 
(Fermionic) electron operator requires multiplying this 
Bose operator by a Jordan- Wigner "string" : 



0(x) 



in > , n(x') 

/—/IE 7 <X V ' 



: j X p(x') 



i(k F x+Q) 



(4.14) 



Since this string operator carries momentum kp, the 
resulting Fermionic operator Oe %Lp should be identified 
with the right moving continuum Fermi field, ip B . We 
have thereby identified the correct Bosonized form for 
the (continuum) electron operators: 



_ e i<t>p(x). 



= <P- 



P0, 



(4.15) 



with P = R/L = ±. From Eqn. [4.10| the chiral Boson 
fields (j>p can be shown to satisfy the so-called Kac-Moody 
commutation realtions: 

\(j> P {x),<j)p(x')] = iPnsgn(x~ x'), (4.16) 

[Mx),M*)] = ™ ■ ( 4 - 17 ) 

These commutation relations can be used to show that 
ip R and ip L anticommute. 

It is instructive to re-express the Bosonized Hamilto- 
nian density in terms of the chiral boson fields, 

H = irv F [n% + n% (4.18) 

where we have defined right and left moving densities 



2tt ■ 



(4.19) 
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which sum to give the total density, nji + til = p- These 
chiral densities can be expressed in terms of the chiral 
electron operators as, 



tppipp := ^ P ipp - (ippipp). 



(4.20) 



Notice that the Bosonized Hamiltonian decouples into 
right and left moving sectors. 

An advantage of Bosonization is the ease with which 
electron interactions can be incorporated. Consider a 
(short-range) density-density interaction added to the 
original lattice Hamiltonian. Using Eqn. i.t this can be 



decomposed into the continuum Dirac fields, and will be 
quartic and spatially local. Due to momentum conserva- 
tion, only three terms are possible: Two chiral terms of 
the form (ippipp) 2 with P = R/L, and a right/left mix- 
ing term of the form, ijr R i}i 



Under Bosonization 



the chiral terms are proportional to (d x <fi P ) 2 , and can be 
seen to simply shift the Fermi velocity in Eqn. 4.18. The 



right/left mixing term also Bosonizes into a quadratic 
form proportional to (d x 9) 2 — (d x (p) 2 . When added to 



the Hamiltonian in Eqn. 4.12, this term can be absorbed 



by shifting both the Fermi velcocity and the dimensionless 
Luttinger parameter, g, which is then no longer equal to 
one. For repulsive interactions g < 1, whereas g > 1 with 
attractive interactions. This innocuous looking shift in 
g has profound effects on the nature of the electron cor- 
relation functions. In fact, it leads to new chiral op- 
erators which have fractional charge, ge. The result- 
ing onepdimensional phase is usually called a "Luttinger 
liquid" B3 For electrons with spin or for Id models with 
multiple bands, the quartic Fermion operators can have 
even more dramatic consequences, for example opening 
up energy gaps as we shall see in Sec. fv|. 

The Lagrangian density in the Bosonized representa- 
tion takes the form of a free scalar field, 



£ = — tti 



(4.21) 



with .9 = 1 for the free Fermion gas, and g ^= 1 in the in- 
teracting Luttinger liquid. The Greek index fj, runs over 
time and the spatial coordinate, /j, = 0, 1 = t, x. Here 
kq = l/irv and ki = —v 2 kq. When re-expressed in terms 
of 9 the Lagrangian takes the identical form, except with 
g —> 1/g for the Luttinger liquid. Changing from the <p 
to the 8 representation can be viewed as a duality trans- 
formation. In Sec. VIII we will consider an analogous 



H = ~tYl [cUx)c a (x r )+h.c.]- (5.1) 

(xx 1 ) x 

where n{x) = c^ a (x)c a (x), and the summation is taken 
over near neighbors on the two- leg ladder, with y = 1,2. 
Due to a parity symmetry under interchange of the two 
legs of the ladder, it is convenient to consider even and 
odd parity bonding and anti-bonding operators: 

b a (x) = -j=[c a (x,y = l) + c a (x,y = 2)], (5.2) 



a a {x) = -j=[c a {x,y = 1) - c a (x,y = 2)], (5.3) 

which depend only on the coordinate x along the lad- 
der. The Hamiltonian splits into even and odd contribu- 
tions, Hq = H a (a) + Hi>(b). Each is a one-dimensional 
tight binding model which can be readily diagonalized by 
transforming to momentum space, 



(5.4) 



and similarly for the anti-bonding operator. Here N de- 
notes the number of sites along the ladder. The diagonal 
form is 



— / J 
k 



k\a\ a a ka + e b k b\ a b ka ], 



(5.5) 



which describes two one-dimensional bands with disper- 
sii 
5. 



sion e^ b = —2t cos k±t — fi. These are sketched in Figure 





1R/ 












~X\1L 


t 2L v^- 





duality transformation in two spatial dimensions. 



Fig. 5: A two-leg ladder and its band structure. In the 
low-energy limit, the energy dispersion is linearized near 
the Fermi points. The two resulting relativistic Dirac 
Fcrmions arc distinguished by pseudospin indices i = 1, 2 
for the anti-bonding and bonding bands, respectively. 



V. 2 LEG HUBBARD LADDER 

A. Bonding and antibonding bands 

We now consider electrons hopping on a two-leg ladder 
as shown in Fig. 5. The kinetic energy takes the form, 



Focussing on the case at half-filling with one electron 
per site (p = 0), both bands intersect the Fermi en- 
ergy, tp = 0. There are four Fermi points at ikFi 
and ±/cf2; f° r the antibonding and bonding bands, re- 
spectively. Gapless particle/hole excitations exist at each 
of the four Fermi points. Due to particle/hole symme- 
try present with near neighbor hopping, e£ + e h k+7r = 0, 
which implies that kp\ + kp2 = tt. Moreover, the Fermi 
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velocity in each band is the same, hereafter denoted as 
v. It is instructive to plot these Fermi points in two- 
dimensional momentum space, taking transverse momen- 
tum k y — 0,7r for the two bands, as shown in Figure 6. 
The four Fermi points can be viewed as constant k y slices 
through a two-dimensional Fermi surface. 

As we shall see, with even weak electron interactions 
present the gapless Fermi points are unstable, and a gap 
opens in the spectrum. Of interest are the properties of 
the resulting Mott insulator. As discussed in Section II], 
for strong interactions mapping to a spin model is possi- 
ble, and the electron spins across the rungs of the ladder 
are effectively locked into singlets: 



\RS) 



1 



V2 



[4(1)4(2) - c|(l) C t(2)]|0), (5.6) 



where y = 1,2 refers to the two legs of the ladder, and 
we have suppressed the rung position x. The state |0) 
denotes a rung with no electrons. 




Fig. 6: Fermi points for the two-leg ladder plotted in 
the two-dimensional Brillouin zone, with the antibonding 
band (denoted 1) at k y = it and the bonding band (2) at 
k y = 0. The shaded region represents the Fermi sea for 
a two-dimensional square lattice model at half-filling. 

It is extremely instructive to re-express this rung- 
singlet state in terms of the bonding and anti-bonding 
operators. One finds, 



\RS) = ^=[b\b\-a\a\]\0), 



(5.7) 



a linear combination of adding a singlet (Cooper) pair 
into the bonding and antibonding orbitals. This paired 
form is suggestive of superconductivity. Indeed, when 
viewed in momentum space, the ground state of a super- 
conductor is a product of singlet pairs with zero center 
of mass momentum at different points around the Fermi 
surface. In an s-wave superconductor, the pairs are all 
added with the same sign, but if the pairs are formed with 
a relative angular momentum (eg. d-wave) sign changes 
are expected. But notice the most important relative mi- 
nus sign in the rung singlet state! The spin-liquid phase 



of the two-leg ladder is evidently related to a paired su- 
perconductor with non-zero angular momentum. Since 
pairing in the bonding band at k y = has a positive sign 
and pairing in the anti-bonding band at k y = 7r is neg- 
ative, in the two-dimensional Brillouin zone (see Figure 
6) the sign is proportional to k% — 
so-called d x 2_ y 2 pairing symmetry. 

If the interactions are weak, it is legitimate to focus 
on electronic states near the Fermi points. As in Section 



ky, consistent with a 



IV, the electron operators can be conveniently decom- 
posed into continuum fields near the Fermi points which 
vary slowly on the scale of the lattice. Denoting cy = a 
and C2 = 6, the bonding and antibonding operators are 
expanded as, 



-ikt. 



(5.8) 



with i — 1,2. Upon linearizing the spectrum around 
the four Fermi points the kinetic energy takes the form, 
Hq = J dxTLo, with Hamiltonian density, 



Ho 



(5.9) 



This Hamiltonian describes massless Dirac Fermions, 
with four flavors labelled by band and spin indices. Im- 
plicit in this theory is a momentum cutoff, A, whose in- 
verse exceeds the lattice spacing. Only modes with mo- 
mentum \k\ < A are included in these continuum fields. 
Since the spectrum is massless, this simple theory is "crit- 
ical" and scale invariant behavior is expected. This can 
be seen by considering the (Euclidian) action, written as 
a space-time integral of the Lagrangian density, 



S = drdxLo, 



(5.10) 



(5.11) 



with P = R/L, and r denoting imaginary time. The 
partition function, Z = Trexp{—(3Ho), can bo-expressed 
as a (coherent state Grassman) path integral,E3 

Z = J [Di)][D$\e- s ^\ (5.12) 

A simple renormalization group can be implementecSllzl 
by first integrating out fields ip(k,uj) with momentum k 
lying in the interval A/6 < |fc| < A, with rescaling param- 
eter b > 1. Since modes with different momentum and 
frequency are not coupled, the action takes the same form 
after this integration, except with a smaller momentum 
cutoff, A/6. The renormalization group transformation 
is completed by a rescaling procedure which returns the 
cutoff to it's original value: 



bx; 



br; V -> 6~ 1/2 i/>. 



(5.13) 



The field rescaling has been chosen to leave the action in- 
variant. This simple theory is at a renormalization group 
fixed point. 



<) 



B. Interactions 



ir i ' T 11 



a i3 1 mj 



(5.18) 



Electron-electron interactions scatter right-moving 
electrons into left-moving electrons and vice-versa. We 
consider general finite-ranged spin-independent interac- 
tions, but assume that the typical interaction strength, u, 
is weak - much smaller than the bandwidth. We focus on 
the effects of the interactions to leading non- vanishing or- 
der in h. In this limit it is legitimate to keep only those 
pieces of the interactions which scatter the low energy 
Dirac Fcrmions. A general four Fermion interaction on 
the two-leg ladder (such as the Hubbard u) can be read- 
ily decomposed in terms of the continuum Dirac fields. 
It is instructive to see how these quartic terms in ip(x) 
transform under the rescaling transformation Eqn. [5.13 



A simple quartic term with no spatial gradients is seen 
to be invariant, so that these operators are "marginal" 
under the renormalization group. The corresponding in- 
teraction strengths will "flow" under the renormalization 
group transformation due to non-linear interaction ef- 
fects. On the other hand, a quartic term involving gra- 
dients such as i/2 (V^ d x ip) 2 , would rapidly scale to zero 
under rescaling: w 2 — * and can thus be ignored. 

Moreover, four-Fcrmion interactions which are chiral, say 
only scattering right movers, do not renorm alizc .to low- 
est order in u and can thus also be neglect edciE3 As dis- 
cussed in Sec. IV, these terms simply lead to small shifts 
in the Fermi velocity. All of the remaining four-Fermion 
interactions can be conveniently expressed in terms of 
currents, defined as 

Iij = ^ia^af}^ 3 p, Iij = ^ ^a( ecr )a/3V'j/3, (5.15) 

where the R, L subscript has been suppressed. Both J 
and I are invariant under global SU(2) spin rotations, 
whereas J and J rotate as SU(2) vectors. Due to Fermi 
statistics, some of the currents are (anti-)symmetric 



(5.16) 



so that la = (no sum on i). 

The full set of marginal momentum-conserving four- 
Fermion interactions can be written 



n {1} , , 



u lj ° RIJ <* LtJ 5 



fijJfii 



J 



(5.17) 



Here and bij denote the forward and backward 
(Cooper) scattering amplitudes, respectively, between 
bands i and j. Summation on i,j = 1,2 is implied. 
To avoid double counting, we set fa = (no sum on 
i). Hermiticity implies &12 = &21 and parity symmetry 
(R <-> L) gives /12 = /21, so that there are generally 
eight independent couplings b^f , b%2 \ ^12^1 an d fi2- 
half-filling with particlc/holc symmetry bn = & 2 2- Addi- 
tional momentum non-conserving Umklapp interactions 
of the form 



are also allowed, (here 1 = 2, 2 = 1). Because the 
currents are (anti-) symmetric, one can always 

choose U12 = U21 for convenience. We also take = 
since la — 0. With particle/hole symmetry there are 
thus just three independent Umklapp vertices, u^, u p l2 , 
and Ui2- Together with the six forward and backward 
vertices, nine independent couplings are required to de- 
scribe the most general set of marginal non-chiral four- 
Fermion interactions for a two-leg ladder with parti- 
cle/hole symmetry at half-filling. 

The renormalization group transformation described 
above can be implemented by L working perturbatively for 
small interaction parametersEjE£l. Upon systematically 
integrating out high-energy modes away from the Fermi 
points and then rescaling the spatial coordinate and 
Fermi fields, a set of renormalization group (RG) trans- 
formations can be derived for the interaction strengths. 
Denoting the nine interaction strengths as gi, and setting 
the rescaling parameter b = 1 + d£ with d£ infinitesimal, 
the leading order differential RG flow equations take the 
general form, 



de.gi — Aijkgjgt, 



(5.19) 



valid up to order g 3 . The matrix of coefficients is 
given explicitly in Ref. [l6| 

These nine coupled non-linear differential equations 
are quite complicated, but can be integrated numerically 
starting with initial values appropriate to a lattice in- 
teraction (such as the Hubbard interaction). This inte- 
gration reveals that some of the couplings remain small, 
while others tend to increase, sometimes after a sign 
change, and then eventually diverge. Quite surprisingly, 
though, the ratios of the growing couplings tend to ap- 
proach fixed constants, which are independent of the ini- 
tial coupling strengths, at least over a wide range in the 
nine dimensional parameter space. These constants can 
be determined by inserting the Ansatz, 



9i{£) 



.9/0 



(id - t) ' 



(5.20) 



into the RG flow equations, to obtain nine algebraic equa- 
tions quadratic in the constants giQ. There are various 
distinct solutions of these algebraic equations, or rays in 
the nine-dimensional space, which correspond to differ- 
ent possible phases. But for generic repulsive interactions 
between the electrons on the two-leg ladder, a numerical 
integration reveals that the floms are essentially always 
attracted to one particular ray.EEl This is the spin-liquid 
phase of interest, which we refer to as a d-Mott phase. In 
the d-Mott phase, two of the nine coupling constants, b^ 
and f^i remain small, while the other seven grow large 
with fixed ratios: 
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&12 = 



fP 
J 12 



2 1 1, 



11 - 2 < 2 = 2 U i2 = 9 > 0. 



(5.21) 
(5.22) 



Once the ratio's are fixed, there is a single remaining 
coupling contant, denoted g, which measures the distance 
from the origin along a very special direction (or "ray") 
in the nine dimensional space of couplings. The RG equa- 
tions reveal that as the flows scale towards strong cou- 
pling, they are attracted to this special direction. If the 
initial bare interaction parameters are sufficiently weak, 
the RG flows have sufficient "time" to renormalize onto 
this special "ray" , before scaling out of the regime of per- 
turbative validity. In this case, the low energy physics, on 
the scale of energy gaps which open in the spectrum, is 
universal, depending only on the properties of the physics 
along this special ray, and independent of the precise val- 
ues of the bare interaction strengths. 



C. Bosonization 

To determine the properties of the resulting d-Mott 
phase, it is extremely helpful to Bosonize the theory. As 
discussed in Sec. |V| the (continuum) electron fields can 
expressed in terms of Boson fields: 



•■pi, 



P0i. 



(5.23) 



with P = R/L = ±. The displacement field 9i a and 
phase field ifi a satisfy the commutation relations 

[ip la {x),6jp{x')] = -iir5ij6 a/ 3Q(x' - x). (5.24) 

Klein factors, satisfying 

{«fa)%} = 2 ¥«i3. ( 5 - 25 ) 

have been introduced so that the Fermionic operators 
in different bands or with different spins anticommutc 
with one another. When the Hamiltonian is Bosonized, 
the Klein factors only enter in the combination, Y = 
Ki|K 1 |K 2 |K2i- Since Y 2 — 1, one can take Y = ±1. Here- 
after, we will put r = 1. 



The Bosonized form for the kinetic energy Eq. pS) 
Ho = J2[(d x 6 ia ) 2 + (d x <p ia ) 2 }, 

ATT ' * 



is 



(5.26) 



which describes density waves propagating in band i and 
with spin a. 

This expression can be conveniently separated into 
charge and spin modes, by defining 



'ip 



9it 



(5.27) 
(5.28) 



and similarly for ip. The y/2 ensures that these new fields 
satisfy the same commutators, Eq. ( 5.24 ). It is also con- 
venient to combine the fields in the two bands into a ± 
combination, by defining 



9^ ± = (9^±9 2 ^)/V2, 



(5.29) 



where fi = p, a, and similarly for ip. 

The Hamiltonian density Ho can now be re-expressed 
in a charge/spin and flavor decoupled form, 



% = f ^[(W,±) 2 + (3^±) 2 ]. 

Z7T * — ' 



(5.30) 



The fields 9 p+ and ip p+ describe the total charge and cur- 
rent fluctuations, since under Bosonization, W Pia ^ Pla 
2d x 9 p+ /n and vP^ pl J) pia = 2d x ip p+ /n. 

While it is possible to Bosonize the interaction Hamil- 
tonians in full generality,E3 we do not reproduce it here. 
In addition to terms quadratic in gradients of the Boson 
fields (as in Ho), the Bosonized interaction consists of 
terms bi-linear in cos 29 and cos2ip. More specifically, 
of the eight non-chiral Boson fields (0 M ± and </? M ± ) only 
five enter as arguments of cosine terms. In the momen- 
tum conserving terms these are 6„±, ip p ^ and ip a -. The 
Umklapp terms also involve the overall charge displace- 
ment field, via cos29 p+ . This can be understood by con- 
sidering how the Boson fields transform under a spatial 
translation, x — > x + Xq- The chiral electron operators 
transform as tppi ~~ > ip P ie lp FiX ° , which is equivalent to 
9ia * 9i a + k Fl x Q . Three of the charge/spin and fla- 
vor fields are thus invariant under spatial translations, 
whereas 9 p+ — > 9 p+ + ttxq. The momentum conserving 
terms are invariant under spatial translations, so cannot 
depend on cos 29 p + . 

The full interacting theory is invariant under spatially 
constant shifts of the remaining three Boson fields - 
ip p +, ifa-+ and #p_. For the first two of these, the conser- 
vation law responsible for this symmetry is readily appar- 
ent. Specifically, the operators exp(iaQ) and exp(iaS' z ), 
with Q the total electric charge and S z the total z- 
component of spin, generate "translations" proportional 
to a in the two fields ip p + and To see this, we note 

that Q — J dxp(x) with p(x) = 2d x 9 p +/ir the momen- 
tum conjugate to ip p +, whereas S z can be expressed as 
an integral of the momentum conjugate to ip a +- Since 
the total charge is conserved, [Q, H] = 0, the full Hamil- 
tonian must therefore be invariant under ip p+ — > ip p+ + a 
for arbitrary constant a, precluding a cosine term for this 
field. Similarly, conservation of S z implies invariance un- 
der ip a+ — * ip a+ + a. 

The five Boson fields entering as arguments of vari- 
ous cosine terms will tend to be pinned at the minima 
of these potentials. Two of these 5 fields, 9 a - and (p a -, 
are dual to one another so that the uncertainty princi- 
ple precludes pinning both fields. Since there are various 
competing terms in the potential seen by these 5 fields, 
minimization for a given set of bare interaction strengths 
is generally complicated. However, along the special ray 
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in the nine dimensional space of interaction parameters 
the nine independent coupling constants can be replaced 
by a single parameter g. The resulting Bosonized theory 
is found to reduce to a very simple and highly symmetri- 
cal form when expressed in terms of a new set of Boson 
fields, defined by 



(0, <p) 2 = (0,(p)a+, 

{9,(p) 4 = (<p,6) p -. 



(5.31) 



The first three are simply the charge/spin and flavor 
fields defined earlier. However, in the fourth pair of fields, 
9 and tp have been interchanged. 

In terms of these new fields, the full interacting Hamil- 
tonian density along the special ray takes an exceedingly 
simple form: TL = TLq + Hi, with 



Z7T * — ' 

n. 



(5.32) 



a 

—4<7 cos 20 o cos 20f,. 



(5.33) 



D. d-Mott Phase 

We now briefly discuss some of the general physical 
properties of the d-Mott phase which follow from this 
Hamiltonian. Ground state properties can be inferred by 
employing semi-classical considerations. Since the fields 
ip a enter quadratically, they can be integrated out when 
the partition function is expressed as a path integral over 
Boson fields. This leaves an effective action in terms of 
the four fields 9 a . Since the single coupling constant g 
is marginally relevant and flowing off to strong coupling, 
these fields will be pinned in the minima of the cosine po- 
tentials. Specifically, there are two sets of semiclassical 
ground states with all 9 a = n a n or all (L_= (n a + l/2)7r, 
where n a are integers. It can be showr£3 that these dif- 
ferent solutions actually correspond to the same physical 
state, so that the ground state is unique. Excitations will 
be separated from the ground state by a finite energy gap, 
since the fields are harmonically confined, and instanton 
excitations connecting different minima are also costly in 
energy. 

Consider first those fields which are pinned by mo- 
mentum conserving interaction terms. Since both O a ± 
fields are pinned, so are the spin-fields in each band, 9 ia 
(i = 1,2). Since d x 9i a is proportional to the z-component 
of spin in band i, a pinning of these fields implies that 
the spin in each band vanishes, and excitations with non- 
zero spin are expected to cost finite energy: the spin gap. 
This can equivalently be interpreted as singlet pairing of 



electron pairs in each band. It is instructive to consider 
the pair field operator in band i: 



A; 



^ Lll =^ Ktl e^^ +e -l (5.34) 



With 9i a w 0, ifip can be interpreted as the phase of the 
pair field in band i. The relative phase of the pair field 
in the two bands follows by considering the product 



a 1 a\ = -r e 



(5.35) 



with r = Ki|Ki|K2|K2i = 1- Since #4 = tp p - the relative 
phase is also pinned by the cosine potential, with a sign 
change in the relative pair field, AjAj < 0, correspond- 
ing to an approximate d-wave symmetry. 

To discuss the physics of the remaining overall charge 
mode {9 P +), it is convenient to first imagine "turning off" 
the Umklapp interactions. After pinning the other three 
fields to the minima of the cosine potentials, the pair field 
operator in band i becomes 



(-iy 



(5.36) 



so that tp p+ is the phase of the pair field. In the absence of 
Umklapp scattering, the Lagrangian for this phase field 
is simply, 



(5.37) 



Being in one-spatial dimension, these gapless phase fluc- 
tuations lead to power law decay of the pair field spatial 
correlation function, A*(x)A(0) ~ l/x n . A true super- 
conductor (for d > 1) exhibits (off-diagonal) long-ranged 
order, and this correlation function would not decay to 
zero even as x — ► 00. But in one-dimension a "supercon- 
ductor" can at best exhibit power law decay, .-since true 
off-diagonal long-ranged order is not possible^ Thus, in 
the absence of Umklapp scattering the 2-leg ladder would 
be a one-dimensional d-wave "superconductor" . 

But what is the effect of the momentum non- 
conserving Umklapp interactions? Once the other three 
fields arc pinned in the minima of the cosine potentials 
in the above Hamiltonian Eqn. 5.33, the Umklapp scat- 
tering terms take the simple form, 



Hu = -12gcos29 p 



(5.38) 



This term tends to pin the field 9 p+ . The pair field phase, 
ip p + , being the conjugate field will fluctuate wildly. These 
quantum flucutations will destroy the power-law Id "su- 
perconducting" phase, leading to an exponentially de- 
caying pair-field correlation function. What is the fate 
of this one-dimensional "quantum disordered d-wave su- 
perconductor" ? 

To see this, one simply has to consider the "dual" rep- 
resentation in terms of the 9 p+ field, rather than ip p +. A 
lattice version of this duality transformation is carried out 
in detail in the Appendix. Alternatively, one can obtain 
the dual theory directly from the Bosonized Hamiltonian 
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Eqn. 5.32. The appropriate Lagrangian dual to Eqn. 5.37 
above, is simply 



which describes gapless density waves. These density flu- 
cutations will be pinned by the Umklapp terms in H U1 
leading to a Mott insulator with a gap to charge excita- 
tions. Since there is also a spin-gap this phase is equiv- 
alent to the spin-liquid, discussed at strong coupling in 
terms of the Heisenberg model in Section III. But we 



now see that this spin-liquid phase exhibits supercon- 
ducitng d-wave pairing correlations, despite being an in- 
sulator. The spin-liquid phase can thus be described as a 
quantum disordered one-dimensional d-wave "supercon- 
ductor" . 

The Euclidian acti on associated with the phase La- 
grangian in Eqn. 5.37 is equivalent to the effective Hamil- 



tonian in the low temperature phase of the classical 2d 
xy model, (with imaginary time playing the role of a sec- 
ond spatial coordinate). The 2d xy model can be disor- 
dered by introducing vortices into the phase of the order 
parameter.!^ Ear this it is convenient to go to a dual 
representation.!^ As shown explicitly in the Appendix, 
the dual represention is equivalent to the P + representa- 
tion, with the strength of the Umkla pp te rm playing the 
role of a vortex fugacity. In Section VIII , we will quan- 
tum disorder a two-dimensional d-wave superconductor, 
and it will be extremely convenient to consider a duality 
transformation - a three dimensional version of the 2d 
9 <-> ip duality discussed here. The resulting nodal liquid 
phase will be particularly simple to analyze in the dual 
representation. 



E. Symmetry and Doping 

Due t o th e hig hly s ymmetric form of the Hamiltonian 
in Eqn. 5.32 and 5.35 , it is possible to make considerable 
further progress in analyzing it's properties. Indeed, as 
shown in Ref. [l6[ under a re-Fermionization procedure 
this Hamiltonian is equivalent to the SO (8) Gross- Neveu 
model,Lj which has been studied extensively by particle 
field theorists. The SO (8) Gross-Neveu model posesses a 
remarkable symmetry known as trialitycj which can be 
used to equate the energies of various excited states. In 
particular, the energy of the lowest excited state with the 
quantum numbers of an electron (charge e and s = 1/2) 
is equal to the energy of the lowest lying spinless charge 
2e exited state (a Cooper pair) . This beautifully demon- 
strates pairing in the insulating d-Mott phase: The en- 
ergy to add two electrons of opposite spin far apart is 
twice as large as the energy to add them into a Cooper 
pair bound state. It turns-aut, moreover, that the Gross- 
Neveu model is integrablecB so it is possible to fully enu- 
merate the energies and quantum numbers of all the low 
energy excited statesES (grouped into SO (8) multisets) 
and compute exactly various correlation functions Ha 



We finally briefly mention the effects of doping the d- 
Mott phase away from half-filling. This can be achieved 
by addi ng a chem ical p otential term to the Hamiltonian 



(5.39) in Eqn. 5J32 and |5.33| , with H p = H - /j,Q, where 



is 



the total electric charge: 



= - a 



(5.40) 



Since the field 6 P + is pinned in the cosine potential by 
the Umklapp interaction terms, H Ul for small [i the den- 
sity will stay fixed at half-filling. Eventually, /j, will pass 
through the Mott charge gap and the density will change. 
This occurs via ir instantons in # p+ , connecting adjacent 
minima of the cosine potential. Each instanton carries 
charge 2e, but no spin, so can be intepreted as a Cooper 
pair. In this doped phase, the Umklapp scattering terms 
will no longer we able to freeze the charge fluctuations, 
and one expects gapless excitations in the density and 
pair field phase, f p +. This doped phase will— exhibit 
power-law d-wave superconducting correlations." 



VI. D-WAVE SUPERCONDUCTIVITY 

We now turn to the case of a two-dimensional super- 
conductor which exhibits a particular type of d-wave 
pairing (denoted d x 2_ y 2) appropriate to the Cuprates. 
Our ultimate goal is to quantum disorder this state to 
obtain a description of the "nodal liquid" . There are two 
main distinctions between the 2d d-wave superconductor 
and it's one-dimensional counterpart considered above. 
Firstly, a 2d superconductor exhibits true (off-diagonal) 
long-ranged order at T = 0. But more importantly, due 
to sign changes in the pair wave function, the d x i_ y i su- 
perconductor exhibits gapless quasiparticle excitations. 
We first briefly review BCS theory which gives one a 
powerful framework to describe d -wav e pairing and the 
gapless quasiparticles. In Section VII below we incorpo- 
rate quantum flucutations of the order parameter phase 
to obtain a complete eff ective low-energy theory of the 
d x - y 2 phase. In Section VII] a dual represention is de- 
rived, and used to quantum disorder the superconductor 
in Section IX. 



A. BCS Theory Re-visited 

It is instructive to briefly review BCS theoryJEI fo- 
cussing on the symmetries of the pair wave function and 
the superconducting order parameter. In particular, it 
is important to emphasize the important distinction be- 
tween the wave function for the center of mass of the 
Cooper pair (often ignored) and the wavefunction for the 
relative coordinate. 

Consider a Hamiltonian expressed as a sum of kinetic 
energy and interaction terms, H = Hq + Hi nt , with Hq 
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given in Eqn. 2.2. We consider a rather general form for 
the electron interactions: 



Hint — — V q (k, fe') C fc+qQ C -fe+q/5 C -fe'+q/3 c fe'+(JQi 

(6.1) 



h.h'q 



which is invariant under global charge U(l) and spin 
SU(2) symmetries. For simplicity Umklapp interaction 
terms have been ignored, so that the crystal momen- 
tum is conserved. The interaction term describes a two 
electron scattering process with 2q the total conserved 
momentum of the pair. For a density-density inter- 
action in real space, such as the Coulomb interaction, 



,(k,k') 



\k — fc'|), so is independent of q. 



Superconductivity within BCS theory requires an at- 
tractive interaction (in the appropriate angular momen- 
tum channel) between electrons. But the bare Coulomb 
interaction is of course strongly repulsive. In traditional 
low temperature superconductors, phonons are believed 
to drive the pairing, inducing a retarded attractive inter- 
action at low energies below the deBye energy. Super- 
conductivity in the high temperature Cuprates is prob- 
ably of electronic origin. In this case, retardation lead- 
ing to an attractive interaction at low energies would be 
due to virtual interactions via high energy electron states 
well away from Ef- These processes-can be studied via 
a renormalization group procedure, £3 which consists of 
"integrating out" high energy electron states, and seeing 
how the remaining interactions between those electrons 
near the Fermi energy are modified. This is precisely 
what we implemented in detail for the two-leg ladder in 
Sec. [v| One thereby arrives at an effective low energy 
theory involving electron states within a small energy 
range of width 2 A around Ep, scattering off one another 
with an effective (or renormalized) interaction potential. 
In the following, we view v q (k,k') as an effective low 
energy interaction. For the two-leg ladder the renormal- 
ized potential is given by putting the nine coupling con- 
tants equal to their values along the special ray. Upon 
Boso nizat ion, the effective potential is given explicitly in 
Eqn. 5.33. More generally, the form of the renormalized 
potential will be constrained by the original symmetries 
of the Hamiltonian. Specifically, time reveral and par- 
ity symmetries imply that v q (k,k') is real, and odd in 
it's arguments: v q (k,k') — V- q (— k, — k'). Hcrmiticity 



implies v q (k, k') 



j(k',k). The summation over mo- 



mentum is now understood to be constrained, involving 
only electron operators with energy in a shell of width 
2 A about E F . 

BCS theory can be implemented by considering the 
operator, 



P^il) = c- k+qa c k+qP} 



(6.2) 



which destroys a pair of electrons, with total momen- 
tum 2q. For k near the Fermi surface, and \q\ << kp, 



[Pk(q), Pk(q')] — for q ^ q' , so that the pair opera- 
tor resembles a boson operator, b(q). By analogy with 
Bose condensation, in the superconducting phase one ex- 
pects a non-zero expectation value for the pair opera- 
tor: (P) 7^ 0. The pair operators entering into Hi nt 
are expressed as P = (P) + 6P, and the fluctuations 
SP = P — (P) are presumed to be small. Upon ignoring 
terms quadratic in SP, Pi-mt can be written (dropping 
additive constants), 



(6.3) 



k.q 



where we have introduced the (complex) superconduct- 
ing order parameter (or "gap"), A, defined as, 



A h (l) = y]v q (k,k'){c_ k > +qa c kl+q p) 



k 



(6.4) 



BCS is a self-consistent mean field theory: The full mean 
field (or quasiparticle) Hamiltonian, H qp = Hq + Hi, 
which depends on A, is employed to compute the expec- 
tation value (c-k> + qa Ck> +qp) ■ Upon insertion in Eqn. 6.4 
one obtains a self-consistent equation which determines 
A - the celebrated BCS gap-equation. Notice that H qp 
is bi-linear in electron operators and hence tractable, al- 
though it does involve "anomalous" terms involving pairs 
of creation or annihilation operators. 

Before carrying through this procedure, it is instruc- 
tive to consider the form for the pair wavefunction which 
follows from a non-zero expectation value of the pair op- 
erator, (P) =/= 0. Consider removing a pair of electrons, 
at positions R± r/2, with R the center of mass position 
and r the relative coordinate. The pair wave function 
can be defined as, 



$ afi (R,r) = (c a (R-r/2) Cf3 (R + r/2)), 



(6.5) 



which depends on the spin of the electrons as well as the 
(center of mass and relative) positions. Upon transform- 
ing the electron operators into momentum space, one 
finds that 



<f> af3 (R,r) 



Q 



,iQRjf,a(3 



$ tt "(Q,r), 



with Q the center of mass momentum and 

^(Q,r) = ^j2^ p (Q/ 2 )y k ' r 

k 



(6.6) 



(6.7) 



Notice that the wavefunction in the relative coordinate, 
involves a Fourier transform with respect to the relative 
pair momentum, k. 

It is also instructive to define a spatially varying su- 
perconducting order parameter by Fourier transforming 
the gap function, Afe(qr): 



Af 0r)=5>^A fc (Q/2). 
Q 



(6.8) 
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In the superconducting phase one can often ignore 
the spatial dependence of the complex order parameter 
Ak(x), and indeed in BCS theory this x dependence is 
dropped. However, if one wishes to include the effects 
of quantum fluctuations (to quantum disorder the super- 
conductor) it is necessary to conside r a s patially varying 
order parameter as discussed in Sec. VII below. 

By analogy with Bose condensation, one expects the 
Cooper pairs to be condensed into a state of zero mo- 
mentum, Q = 0. This requires 



( P k (?)) = &q.o{ c ~kaCk(}) 



(6.9) 



which gives a relative pair wavefunction, $(r) = Q(Q = 
0, r) of the form, 

*°'V) = ^£ e * r *Z / '; = <C-*a<*0>. (6.10) 



Due to the electron anticommutation relations one has 



a/3 



<3?_k, which implies that the pair wavefunction 



is antisymmetric under exchange of the two electrons: 
$ a/3 (r) = -<J>^ a (-r). 

When the Cooper pairs are condensed into a state with 
zero momentum, the superconducting order parameter 



becomes spatially uniform: A?^(x) = A™/ 3 , as seen from 



.a/3 



Eqn. 5.4. The mean field Hamiltonian then takes a rather 



simpler form: 



Hi 



9 £l C fco 



-k0 L 



+ h.c], 



(6.11) 



whereas the self-consistentcy condition becomes, 

= - £ ««,(*, fe')<C-fe'aCfe'/3). (6.12) 



Since the full model has a conserved SU{2) spin sym- 
metry, the relative pair wavefunction can be expressed 
as the product of an orbital and a spin wavefunction: 
Q^f = Q /3$fe. The spin piece can be chosen as an eigen- 
function of the total spin of the pair, that is a singlet 
with 5 = or a triplet with 5=1. In conventional 
low temperature superconductors and in the Cuprates 
the Cooper pairs are singlets with, 



4>af3 = $ai$/3] ~ S^Spi, 



(6.13) 



in which case the orbital wavefunction is symmetric: 
= ff-fc = (c_fcjCfcj). (In the superfluid phases of 
3 — He on the other hand, the Cooper pairs have 5 = 1.) 
The superconducting order parameter is then also a sin- 
glet; A^ = 4> a pAk, with A fc = A_ fc satisfying 



Afc = "^£ w o(fc,fc')(c_fc'|Cfc'1 



(6.14) 



For singlet pairing, the final mean field (quasiparticle) 
Hamiltonian becomes, H qp = Hq + Hi with, 



H x =£[A fc c 



fcT -fei 



A fe c -fcl c fcT]- 



(6.15) 



To complete the self-consistency requires diagonalizing 
the quasiparticle Hamiltonian. This is usually doue.in a 
way which masks the spin rotational invarianceO We 
prefer to keep the spin rotational invariance explicit, by 
defining a new set of Fermion operators, for k y > 0: 

Xla(fe) = C ka ; X2a{k) = i<TapC'_ h p, (6.16) 

which satisfy canonical Fermion anti-commutation rela- 
tions: 



(6.17) 



The first index a,b — 1,2 acts in the particle/hole sub- 
space. The o~ y in the definition of xia has been intro- 
duced so that these new operators transform like SU (2) 
spinors under spin rotations: \ aa — > U a pXa/3, with 
U = cxp(i# • er) a global spin rotation. 

In these variables, the quasiparticle Hamiltonian be- 
comes 

H 9P =Y,'xHk)[T g € k + T + A k + T-Al]x(k), (6.18) 
k 

where the prime on the summation denotes over k y pos- 
itive, only, and we have introduced a vector of Pauli ma- 
trices, T a b acting in the particle/hole subspace. Also, we 
are employing the notation r = (t x ± ir v )/ 2. To evalu- 
ate the self-consistency condition Eqn. |6.14| we need the 
anomalous average of two electron fields (the orbital piece 
of the relative pair wavefunction), which is re-expressed 
as, 

\ E( c Tfclc ±feT ) = l -{xHk)r + X {k)). (6.19) 

Diagonalization is now achieved by performing an 
SU(2) rotation in the particle/hole subspace, by defin- 
ing rotated Fermion fields: x(fc) = U(k)x(k), with 



U(k) 



T . Assuming for simplicity that A k is real, 



the appropriate rotation is around the y-axis by an angle 
6 k , U(k) = e~' WkTV / 2 , with 



sin(0 k ) = — ; 

&k 



Ek — 



(6.20) 



In terms of the rotated Fermion fields, x, the quasiparti- 
cle Hamiltonian is diagonal, 



H„ 



Ek X \k)T z x(k) 



k 



(6.21) 



with Ek the quasiparticle energy. Finally, we define a set 
of rotated electron operators via 

Xia(k) = a ka ; X2a(k) = icr^aL^, (6.22) 
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and the quasiparticle Hamiltonian can be re-expressed in 
standard form, 



k 



Ek a ka aka > 



(6.23) 



where we have dropped an additive constant. Notice that 
the quasiparticle energy E k > for all momentum. On 
the Fermi surface, t k = and the quasiparticle energy is 
given by |Afe| - the energy gap. 

To complete the self-consistentcy, the anomalous elec- 
tron avera ge (or relative orbital pair wavefunction from 
Eqn. |6.19| ) is expressed in terms of the quasiparticle 
operators. Upon using the fact that [£^T+[/] diQg 
sin(6)r z /2 one obtains, 



A 



2Ek 



a ka a ka ) 



(6.24) 



which reduces to = —Ak/2Ek at zero tempera- 
ture. At finite temperature the number of quasiparticles 
is simply a Fermi function: (a^aka) = 2f(Eh), with 
f(E) = [exp(f3E) + l] -1 . One thereby obtains the cele- 
brated BCS gap equation: 



A * = -\ E v °( k > k '^ [1 2f{Ek ' )] 



(6.25) 



Recent experimentEEB has established that in the high 
temperature superconductors the orbital pairing symme- 
try is a particular form of d-wave, usually denoted as 



(1^2 — ? 



Here x and y refer to the directions along the 



crystalline axis of a single Cu-0 sheet, within which the 
Cu atoms form a square lattice. In terms of the corre- 
sponding two dimensional momentum, k = (k x ,k y ), the 
angular dependence of the gap function in this state is 
Afc ~ k^ — ky, and from Eqn. 6.24 the orbital piece of the 
relative pair wave function has the same d-wave symme- 
try. 

A novel feature of the d x i_ y i state is that the gap func- 



_ tion Afe vanishes along lines in k — space with k x = ±k y 
corresponding to nodes in the relative pair wave func- 
tion. These lines intersect the (two-dimensional) Fermi 
surface at four points in momentum space. Near these 
four points (or "nodes") in momentum space there are 
electronic excitations with arbitrary low energy, in strik- 
ing constast to the fully gapped s-wave case. These low 
energy quasiparticle excitations dominate the physics of 
the d x 2_ y 2 superconductor at temperatures well below 
T c , leading to power law temperature corrections in such 
quantities as the electronic specific heat and the magnetic 
penetration length. 



C. Continuum description of gapless quasiparticles 



B. d-wave Symmetry 

In a system with rotational invariance the orbital 
piece of th e pair wavefunction, proportional to A& from 
Eqn. 6.24, can be chosen as an eigenstate of angular mo- 
mentum, a spherical harmonic YJ m in three dimensions. 
The simplest case is s-wave, with A& a constant over 
the (spherical) Fermi surface. Real materials of course 
do not share the full continuous rotational symmetry of 
free space. Nevertheless, a superconductor in which A& is 
everywhere positive over the Fermi surface is (loosely) re- 
ferred to as having s-wave pairing - a property of all con- 
ventional low temperature superconductors. Since |Afc| 
is the quasiparticle energy on the Fermi surface, there 
are no low energy electronic excitations in an s-wave su- 
perconductor - the Fermi surface is fully gapped. Within 
BCS theory the magnitude of the (zero temperature) en- 
ergy gap is related to the superconducting transition tem- 
perature: 2|A| ss 3.5fcsT c . The presence of an energy gap 
leads to thermally activated behavior for various low tem- 
perature properties, such as the electronic specific heat 
and the magnetic penetration length. 

It is clear from the self-consistent gap equation 
Eqn. |6.25 that a purely repulsive effective interaction, 
v (k 7 k') > 0, precludes s-wave pairing within BCS the- 
ory (since 1 — 2f(Ek) > 0). In conventional supercoiu 
ductors, phonons are believed to drive s-wave pairing,Lj 
generating an effective attractive interaction at low ener- 
gies. 



It is convenient to obtain a continuum description of 
the gapless d-wave quasiparticles, analogous to the Dirac 
theory desription of the low energy properties of the Id 
free Fermions employed in Section |l^. A continuum 
form can be obtained directly from the general quasi- 
particle Hamiltonian Eqn. 6.18 by specializing to d x i_ y i 



symmetry and then focussing on those momenta close to 
the four nodes where the quasiparticle energy E k = 
(see Figure 7). For a model with particle/hole symmetry 
£fe = — efe+Tr, which together with parity symmetry im- 
plies that the four nodes occur at the special wavevectors 
±Kj, with K x = (7r/2,7r/2) and K 2 = (-tt/2, tt/2). It 
is convenient to introduce two continuum fields tyj, one 
for each pair of nodes, expanded around ±K\, ±K 2 '- 



(6.26) 



Here, the wavevectors q are assumed to be small, within a 
circle of radius A around the origin. With this definition, 
the particle/hole transformation is extremely simple, 



(6.27) 



For this reason it is convenient to always define the con- 
tinuum fields i& around ±Kj , and account for deviations 
of the node momenta from these values by a particle/hole 
symmetry-breaking parameter A. 
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Fig. 7: In the d x 2_ y 2 superconductor the quasiparticle 
energy vanishes at four points (±K\ and i-K^) in the 
Brillouin zone. The dotted line represents the Fermi sur- 
face. The wavevector q is rotated with respect to the a 
and b axis of the square lattice. 

Once we have restricted attention to the momenta near 
the nodes, it is legitimate to linearize in the quasiparti- 
cle Hamiltonian. The resulting theory is more conve- 
niently written in coordinates perpendicular and paral- 
lel to the Fermi surface, so we perform the rotation via 
x — > (x — y)j\pl and y — > (x + y)/y/2, correspondingly 
transforming the momenta q x and q y (see Figure 7) . Lin- 
earizing near the nodes, we put ejti+q = v F qx where v F 
is the Fermi velocity and 



E 1 {q) = \J (vFq x + A) 2 + |A| 2 g 2 , 



(6.32) 



and a similar expression is obtained for Ei except with q x 
and q y interchanged. As usual in Dirac theory, the nega- 
tive energy single particle states with energy —Ej(q) are 
filled but positive energy holes states can be created. As 
expected, the quasiparticle energy vanishes at the nodes 
(q = with particle/hole symmetry A = 0), so the "rel- 
ativistic" particle is massless. Notice that non-zero A 
indeed shifts the positions of the nodes. 

In this continuum description A serves as a complex 
superconducting order parameter for the d x 2_ y 2 state. 
Indeed, when A = one recovers the metallic Fermi sur- 
face and the quasiparticle Hamiltonian describes gapless 
excitations for all q y . Below we will include quantum 
fluctuations by allowing A to depend on space and time. 
Before doing so, it is convenient to see how A transforms 
under a particle/hole transformation. From the transfor- 
mation properties of the electron fields one deduces that 
the gap transforms as, A& — * — Al fc+7r , which is equiv- 
alent to complex conjugation for the (linearized) order 
parameter, 



A -> A*. 



(6.33) 



Together with Eqn. 6.27 this implies that the quasiparti- 
cle Hamiltonian in Eqn. 
metric: — > Hq,. 



6.3C is indeed particle/hole sym- 



A 



Kt+q 



Aq y + 0{q 2 



(6.28) 



where A has dimensions of a velocity. An identical lin- 
earization is possible around the second pair of nodes, 
except with q x <-> q y . It is finally convenient to Fourier 
transform back into real space by defining, 



* 1 



(6.29) 



where the momentum summation is for q < A. The 
continuum fields *£(x) describe long lengthscale varia- 
tions of the quasiparticles, on scales greater than A -1 . 
We thereby arrive at a compact form for the continuum 
quasiparticle Hamiltonian in a d x 2_ y 2 superconductor: 
H qp — + H\ with 

= v\[v F T*id x + (Ar+ + AV _ )ia„]*i 
+ (1 <-> 2; x <-> y), 



(6.30) 



and the particle/hole symmetry breaking term, 



A*t T z ^,-. 

J J 



(6.31) 



The quasiparticle Hamiltonian takes the form of (four) 
Dirac equations in 2 + 1 space-time dimensions, and can 
be readily diagonalized. For the first pair of nodes one 
obtains the relativistic dispersion, 



VII. EFFECTIVE FIELD THEORY 

A. Quasiparticles and Phase Flucutations 

Our goal in this section is to obtain a complete low- 
energy effective theory for the d x 2_ y 2 superconductor. 
This task is complicated by the existence of additional 
gapless excitations, besides the quasiparticles. Specifi- 
cally, since the global U(l) charge conservation symmetry 
(c a — > e l6o c a ) is spontaneously broken by the existence 
of a non-zero order parameter, A/0, gapless Goldstone 
modes are expected. (In a three-dimensional supercon- 
ductor these modes are actually gapped, due to the pres- 
ence of long-ranged Coulomb interactions, but would be 
gapless for a thin 2d film.) These modes propogate in the 
phase of the complex order parameter. Thus a correct low 
energy theory for the d x 2_ y 2 state requires consideration 
of a spatially varying order parameter, A(x). Generally, 
both the magnitude and the phase of the complex order 
parameter can vary, but we will focus exclusively on the 
phase fluctuations, writing 



A{x) = \A\e 



(7.1) 



with |A| a (real) constant. Since amplitude fluctuations 
are costly in energy this should suffice in the supercon- 
ducting phase, and will also allow us to describe the nodal 
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liquid in which superconductivity is destroyed by phase 
fluctuations. The desired low energy effective theory can 
be obtained from symmetry considerations, and the form 
of the continuum quasiparticle Hamiltonian. A more mi- 
croscopic approach, discussed briefly below, would entail 
integrating out high energy degrees of freedom in a func- 
tional integral representation. 



1. Symmetry considerations 

Since the BCS gap equation has a degenerate mani- 
fold of solutions for arbitrary phase ip, the energy should 
only depend on gradients of ip(x). The appropriate La- 
grangian which describes the fluctuations of the phase of 
the superconducing order parameter can thus be devel- 
oped as a gradient expansion, with lowest order terms of 
the form, 



where the Greek index /z runs over time and two spatial 
coordinates: u — 0,1,2 = t,x,y. Here kq is equal to 
the compressibility of the condensate (ignoring for the 
moment long-ranged Coulomb forces) and Kj — V^, Kq 
(for j = 1,2 = x,y) with v c the superfluid sound veloc- 
ity. This form is largely dictated by symmetry. Parity 
and four-fold rotational symmetry determine the form of 
the spatial gradient terms. The stiffness coefficients, k^, 
can be estimated as follows. The pair compressibility kq 
should be roughly one half the electron compressibility - 
at least for weak interactions. If the pairing is electronic 
in origin, the Fermi velocity sets the scale for v c . 

In general a Berry's phase termu linear in dt<p is al- 
lowed, 



jC-Berry — ngOtip, 



(7.3) 



where no is a two-dimensional number density. For a 
model with particle/hole symmetry which must be in- 
variant under 



-f, 



(7.4) 



(which follows from the particle/hole transformation 
properties of the order parameter A ~ e %v — ► A*) it 
naively appears that the nuipber density no must vanish. 
However, this is not the caseO-3. To see this it is necessary 
to return to the lattice where the term in the (Euclidian) 
action which follows from Cserry is, 



S 



Berry ^?7o 



liTi 



dr^^dripi 



(7.5) 



where i labels the sites of a square lattice with lattice 
spacing set to one and (3 — 1/fcgT . The partition func- 
tion is expressed as a functional integral of exp(-S) over 



configurations ¥?i(t), with /3 periodic boundary condi- 
tions on the fields e %{f '. This implies the boundary condi- 
tions, 



<p i (0) = <Pi(O)+2irN i 



(7.6) 



with integer winding numbers Nj . We thus see that the 
Berry's phase term contributes a multiplicative piece to 
the partition function (in each winding sector) of the 
form; 



exp( Snerry) 



_ i2nn N w 



(7.7) 



with integer Nw = XI - Nj. Under the particle/hole 



transformation Eqn. 7.4, the winding numbers change 
sign, Nw — * — Nw- The Berry's phase term is thus in- 
variant under the particle/hole transformation provided 
no is integer or half-integer. 

The appropriate value for no can be readily determined 
by obtaining the lattice Hamiltonian associated with the 
(7 2) Lagrangian density C v + L Berry The first contribution 
can be conveniently regularized on the lattice as, 

L v = -t ^ cos((pi - (fj) - - ^(dtfi) 2 . (7.8) 

<i,j> i 

Upon inclusion of the Berry's phase term this gives the 
lattice Hamiltonian, 

H v = —t ^2 C0S {fi - fj) + U X^" 1 _ U °) 2 - 

<i,j> i 

Here rn denotes a Cooper-pair number operator, canoni- 
cally conjugate to the phase fields: 



[fi 



iSi 



(7.10) 



The first term in H v describes the hopping of charge 
2e (spinless) Cooper pairs between neighboring sites of 
the lattice, and the second term is an onsite repulsive 
interaction. The parameter no plays the role of an "off- 
set" charge and determines the average number of Cooper 
pairs per site. For the Hubbard model at half-filling with 
one electron per site, the number of Cooper pairs clearly 
equals one-half the number of lattice sites. This is espe- 
cially apparent in the limit of very large attractive Hub- 
bard interaction when the electrons pair into on-site sin- 
glets, but is expected to be more generally valid. Thus, 
it is clear that one should take no = 1/2. Tuning away 
from particle/hole symmetry with a chemical potantial 
11, corresponds to changing no away from one-half. 

In the superconducting phase one expects that the 
winding numbers will all vanish, since the phase field (p 
is essentially constant in space and time, and the Berry's 
term plays no role. But when the superconductor is 
"quantum disordered" , the phase field fluctuates wildly 
with signifigant winding, and inclusion of the Berry's 
phas e term is expected to be important (but see Section 
VLB] below). 
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It remains to couple these phase fluctuations to the 
gapless quasiparticles. Since the order parame ter A di- 
rectly enters the quasiparticle Hamiltonian Eqn. |6.30| , one 
can readily guess the appropriate coupling. We should 
simply replace A — > v A e 4V with w A real. Since ip varies 
spatially, some care is required. In the quasiparticle 
Hamiltonian we let, 



^' 2 {id y )e^' 2 , 



(7.11) 



and similarly for the t~ term. This "symmetric" form 
leads to an hermitian Hamiltonian, physical currents, and 
respects the symmetries of the problem. A careful deriva- 
tion of Eq. [7.1l| is given below. With this prescription, 
the quasiparticle Hamiltonian becomes 



H 



qp 



E ,T,t 

s=± 



*\[vrr*id x + v & T s e tsv/2 (td y )e tslp/2 }^i 



+ (1 «-> 2; x «-> y). 



(7.12) 



Since tp can also fluctuate with time, it will convenient to 
consider the time dependence via a Lagrangian formula- 
tion. The Lagrangian density is 



: qp = v]idt9j - n qp . 



(7.13) 



The full low-energy effective Lagrangian in the d-wave 
superconductor is obtained by adding the two contribu- 
tions: C v + C qp . 



2. Microscopic approach 

To illustrate how one might try to "derive" this effec- 
tive theory from a more microscopic starting point, we 
briefly consider a simple model Hamiltonian, 



H = Ho — V 



cl l (x)cl(x')c (x')c Q (x), (7.14) 



where Hq is the usual kinetic energy describing hopping 
on a 2d square lattice and we have added an attractive 
near-neighbor interaction with strength V. To derive the 
effective field theory, it is convenient to express the parti- 
tion function Z = Tx£~" H , as an imaginary time coherent 
state path integral,cj 



Z 



Vc][Vc}e~ 



(7.15) 



where c and c are Grassman fields and the Euclidean 
action is simply 



Z = J[Vc\[Dc}[DA][VA*}e- s \ (7.17) 

with S± = JdT[J2gC a (x)d T c a (x) + H e fj]. The effective 
Hamiltonian can be decomposed into H e tf = Ho+Hi nt + 
H A , with 



Hint — 



{xx') 



[A%c a {x)cp{x>) + 



h.c. 



- y \A a J,\ 2 . 

y / j I xx' I 

(xx') 



(7.18) 



(7.19) 



Eqs. 7. 18 - 7. IE form a basis for studying the original 
electron model. At this stage BCS mean field theory 
could be implemented by integrating out the electron de- 
grees of freedom to obtain an effective action only de- 
pending on A, S e ff(A). Minimizing this action with 
respect to A would give the gap equation. One could 
imagine including fluctuations by expanding about the 
saddle point solution. But for a d-wave superconduc- 
tor this procedure is problematic, since integrating out 
gapless quasiparticles will generate singular long-ranged 
interactions in S e ff(A). It is preferable to retain the 
gapless quasiparticles in the effective theory, and only 
integrate out the high frequency electron modes which 
are well away from the nodes. In this way, the dynam- 
ics and interactions generated for the order parameter A 
will be local. 

Rather than trying to implement this procedure, we 
content ourselves with arguing that the "symmetric" pre- 
scription adopted above indeed gives the correct form for 
the phase-quasiparticle coupling term. To this end we 
focus on singlet pairing, defining 



A^=A(x,f)(^ T 5 m 



(7.20) 



The triplet pieces of A are presumed to be massive, so 
that they can be safely integrated out. Since A lives on 
the bonds, it is convenient to associate two such fields 
with each site on the square lattice, i.e. 



Ai(f) ee A(x,x + ei), 
A 2 (x) ee A(x,x + e 2 ), 



(7.21) 
(7.22) 



where e\, £2 are unit vectors along the a and b axes of the 
square lattice, respectively. The interaction Hamiltonian 
becomes, 



S= j dr \ ^c a (x)d T c a (x) + H[c,c] 



(7.16) 



We consider here only T = 0, for which the r integration 
domain is infinite. The possibility of a d-wave supercon- 
ducting phase can be entertained by decoupling the above 
action using a Hubbard-Stratonovich transformation: 



Hint — 



Aj(x) c|(a?)cj(x + e.j)— t-*-*! 



(7.23) 



3i x 



where the sum includes all lattice sites and j = 1,2. The 
magnitudes of Ai and A2, as well as their relative sign, 
are determined by the effective action generated upon 
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integrating out the high-energy modes. For a d-wave su- 
perconductor the effective action will be minimized for 
Ai = — A2 = Aoe l(p , up to massive modes. We can now 
take the continuum limit. For agreement with Sec. [v| , 
we define w A = 2V2A , or Ai = -A 2 = A/2\/2. In ad- 
dition, we take the continuum limit of the electron fields, 
using the decompositions 



+ ^l 2 (-0^ + *221^, 



and the hermitian conjugates of these equations. In- 
serting these into Eq. [7.23 , gradient-expanding the W 
fields, and rotating 45 degrees to x — y coordinates along 
the (71", 7r) and (— 7r,7r) directions, one obtains Hint — 
J (PxHint, with 



'Hint — 



h.C. 



+ (1 <-> 2, x «-> y). 



(7.24) 



This form is identical to the A term in Eq. 6.30| when 
the order parameter A is constant, but the symmetric 
placement of derivatives is important in the presence of 
phase gradients. In particular, now let A = v & e l(p and 
integrate by parts to transfer the derivative in the second 
term from the W to the A\l/ combination. Upon using 
the operator identity 

i (e lip id y + idye^) = eWidyeW, (7.25) 

this becomes identical to the symmetrized form o f the 
phase-quasiparticle interaction hypothesized in Eq. 7.12 



B. Nodons 



transformation in Eqn. 2.8). As we shall see, when the 



d-wave superconductivity is quantum disordered, these 
new fields will play a fundamental role, describing low 
energy gapless excitations, centered at the former nodes. 
For this reason, we refer to these fermions as nodons. For 
completeness, we quote the symmetry properties of the 
nodon field under a particle/hole transformation. Since 
ip — » — ip, one has simply 



i\) — ► ip^ . 



(7.27) 



The full Lagrangian in the d-wave superconductor, 
C = C ip + C qp , can be conveniently re-expressed in terms 
of these nodon fields since C qp = + Ci nt + C\ with a 
free nodon piece, 

C^f, = ip\[id t - v F r z id x - v^idyjip! 

+(l<-> 2,a: <->!/), (7.28) 

interacting with the phase of the order-parameter: 

Cint = duipJf,. (7.29) 

Here the electrical 3-current J ^ is given by 



vp 1 1 / 



(7.30) 



(7.31) 



Because the transformation in Eq. [7.26 is local, identical 
expressions hold for these currents in terms of the quasi- 
particle fields, ty. The form of the particle/hole asymme- 
try term remains the same in terms of the nodon fields: 



(7.32) 



Treatment of quantum phase fluctuations is compli- 
cated by the coupling between the quasiparticle Fermion 
operators, '5, and expon entia ls of the phase tp, as seen 
explicitly in TL qp in Eqn. 7.12, The form of the coupling 
is determined by the electric charge carried by which 
is uncertain - being built from electron and hole opera- 
tors. To isolate the uncertain charge of ^ it oh extremely 
convenient to perform a change of variablesjlj defining a 
new set of fermion fields ipj via 



(7.26) 



In the superconducting phase, and in the absence of 
quantum flucutations of the order-paramater phase, one 
can set (p = 0, and these new fermions are simply the 
d-wave quasiparticles. However, when the field ip is dy- 
namical and fluctuates strongly this change of variables 
is non-trivial. In particular, the new fermion fields ip are 
electrically neutral, invariant under a global U(l) charge 
transformation (since <p — * p + 29o under the U(l) charge 



It is instructive to re-express the components of the 
currents J u back in terms of the original electron opera- 
tors. One finds 



Jn = i(c t 



(7.33) 



(with an implicit spin summation) which corresponds 
physically to the total electron density living at the 
nodes, in units of the Cooper pair charge. Similarly, 



J •-— fr 1 c 



K 



) 



(7.34) 



corresponds to the current carried by the electrons at 
the nodes. Thus, J u can be correctly interpreted as the 
quasiparticles three-current. 

To complete the description of a quantum mechani- 
cally fluctuating order parameter phase interacting with 
the gapless fermionic excitations at the nodes, we min- 
imally couple to an external electromagnetic field, A^. 
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Since the nodon fermions are neutral, the only cou- 
pling is to the order-parameter phase, via the substi- 
tution d a p — > d^ip — 2A^. Here we have set the elec- 
tron charge e = 1, with a factor of 2 appropriate for 
Cooper pairs. The final Lagrangian then takes the form 



(7.40) 



£ £ip £~*ih ~\~ £jr. 



C\, with 



Hi 



(d^ - 2A ll )J li 



(7.35) 



(7.36) 



and Cif, still given by Eq. 7.28. Here we have dropped 



the Berry's phase term, which is not expected to play 
an important role in the superconducting phase. Long- 
ranged Coulomb interactions could be readily incorpo- 
rated at this stage by treating Ag as a dynamical field 
and adding a term to the Lagrangian of the form, C cou i = 
(\/2){djA ) 2 . The spatial components of the electromag- 
netic field, Aj, have been included to keep track of the 
current operator. 



1. Symmetries and Conservation Laws 

If the full effective Lagrangian £ is to correctly describe 
the low energy physics it must exhibit the same symme- 
tries as the original electron Hamiltonian - the most im- 
portant being charge and spin conservation. Since the 
tp operators are electrically neutral the full U(l) charge 
transformation is implemented by ip — > ip + 20$ for con- 
stant #o> an d £ is indeed invariant. Moreover, the La- 
grangian is invariant under tp a — > U a pil>/3 for arbitrary 
(global) SU(2) spin rotations U = exp(i0 ■ a). Since the 
Cooper pairs are in spin singlets, all of the spin is carried 
by the nodons. 

As usual, associated with each continuous symmetry 
is a conserved "charge" which satisfies a continuity equa- 
tion (Noether's theorem). Since the Lagrangian only de- 
pends on gradients of <p, the Euler-Lagrange equation of 
motion reduces to the continuity equation, 



W - 0, 



(7.37) 



where the total electric 3-current is given by J 1 " 1 = 
djC/d(d„ip) = -dC/dA^. This gives, 

4 ot = Kfl (d^ - A^) + (7.38) 

where the first term is the Cooper pair 3-current and the 
second the quasiparticles current. 

The analogous conserved spin currents can be obtained 
by considering infinitesimal spin rotations, 



U = l+id(x,t) ■ cr, 



(7.39) 



for slowly varying 0(x,i). Under this spin rotation the 
Lagrangian transforms as, 



with given below. After an integration by parts, in- 
variance of the action S under global spin rotations im- 
plies continuity equations duju — for each of the three 
spin polarizations, j . The space-time components of the 
conserved spin currents are given explicitly by, 



Jo 



+ (1 



2), 



(7.41) 



(7.42) 



and j y the same as j x except with tp± <-> tp 2 - Notice 
that in contrast to the electrical current, the spin cur- 
rent operator has a contribution which is proportional to 
the velocity tangential to the Fermi surface, d a , which 
is anomalous when re-expressed in terms of the original 
electron operators. 

Surprisingly, the effective Lagrangian exhibits addi- 
tional continuous symmetries, not present in the origi- 
nal Hamiltonian. Firstly, £ is invariant under separate 
SU(2) spin rotations on the two pairs of nodes, tpj for 
j = 1,2. Moreover, the Lagrangian is also invariant un- 
der two additional U(l) transformations tpj — ► e tej tpj for 
arbitrary constant phases, 8j. These latter symmetries 
imply two new conserved "charges", tpj tpj (no sum on 
j). We refer to these conserved quantities as "nodon 
charges". The associated conserved nodon 3-currents 
take the same form as the spin currents above, except 
replacing a/2 by the identity. As seen from Eq. [7.31 , the 
conserved nodon charges are proportional to the quasi- 
particle electrical current, since Jj = (vp/2)tpjtpj. 

It is possible to add to C additional interaction terms 
which are consistent with the original U(l) and SU(2) 
symmetries, but do not conserve the "nodon charge" . 
Specifically, anomalous quartic interaction terms of the 
form tp 4 arise from Umklapp scattering processes in the 
original electron Hamiltonian and clearly change the 
nodon charge. However, such interactions are unimpor- 
tant at low energies due to severe phase space restric- 
tions. To see this, consider how the action, S = J d 2 xdtC 
transforms under a renormalization group (RG) rescaling 
transformation, 



bx„ 



tp -> 



LP 



-1/2 



p. 



(7.43) 



with rescaling parameter b > 1. By construction, this 
leaves the quadratic pieces and S v invariant, but in- 
teraction terms such as utp 4 scale to zero under the RG 
(b — > oo) since u — > u/b. It is the T = "fixed point" the- 
ory described by the quadratic terms which exhibits the 
additional symmetries. Incidentally, the coupling term 
Li n t above also scales to zero (as 6 -1 / 2 ) under the renor- 
malization group. In the resulting quadratic theory the 
quasiparticles and phase fluctuations actually decouple. 
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2. Superfluid stiffness 

The above effective theory is particularly convenient 
for examining very low temperatures properties of the 
d x 2_ y 2 state. Of interest are charge response functions 
such as the electrical conductivity and the superfluid stiff- 
ness (measureable via the penetration length). The spin 
excitations (carried by the quasiparticles) can also be 
probed via resonance techniques, such as NMR and ESR. 
Impurity scattering can be readily incorporated by cou- 
pling a random potential to the electron density (which 
can be re-expressed as a nodon bi-linear). 

For illustrative purposes we briefly consider the quasi- 
particle contribution to the low temperature superfluid 
stiffness and extract the famous T-linear dependence. 
For a Galilean invariant system of mass m bosons the 
superfluid stiffness K s equals the superfluid density di- 
vided by m. But more generally K s can be extracted 
rather directly by considering the. response of the system 
to a transverse vector potential.ES We set Aq — and de- 
compose the static vector potential Aj into longitudinal 
and transverse pieces: 



Aj = Afj + A 



(7.44) 



with djAt.j = and e-yA^j = 0. The superfluid stiffness 
is then given by, 



1 d 2 F 
VdA% : 



(7.45) 



where F = —ksTlnZ is the Free energy and V — > oo 
is the area of the 2d system. Here At >x can be taken 
spatially constant. 

To extract F the partition function can be written as 
an imaginary time coherent state path integral,c3 



Z= [Vtp][V^][V^]exp{-S B ), 



(7.46) 



with Euclidian action S E = J d 2 xdrC E - The longitudi- 
nal vector potential, which can be expressed as a gradi- 
ent of a scalar held Aij = djA, can be eliminated en- 
tirely by shifting tp — > tp + A. Moreover, the crossterm 
between djtp and A t j vanishes since A t is divergence- 
less. The Gaussian integral over tp can then be read- 
ily perfomed and simply generates an irrelevant inter- 
action term (J ~ (ip^ip) 2 ) which can be ignored. One 
thereby arrives at an effective action depending only on ip 
and Aj with associated Hamiltonian density of the form: 
Ti-eff = 7~Lip + Ha, with Ti^ the free nodon Hamiltonian 
and 



±K s Al j+ A t jJj. 



(7.47) 



Here K° = kqV 2 is the superfluid stiffness from the 
Cooper pairs, and Jj = (v F /2)ip^tpj. Notice that the 



(transverse) vector potential acts as an effective chemi- 
cal potential for the "nodon charge" density, p n = ipjipj- 
Thus, the superfluid stiffness can be expressed in terms 
of the nodon "compressibility" as 



K° s - (v F /2) 2 Kn , 



(7.48) 



where K n = dp n /dfj, n and \i n — (v F /2)A t . x is the nodon 
"chemical potential" . 

The nodon compressibility can be extracted by diag- 
onalizing the Hamiltonian, H^,. From the first pair of 
nodes one obtains the free Fcrmion form, 



(7.49) 



where E\{q) is given in Eqn.|6~32|and we have suppressed 
the spin index. Here a and b are particle and hole oper- 
ators, respectively. The nodon charge is simply, 



(7.50) 



where the averages are taken with TL^ — fi n Pn- At finite 
temperatures one obtains 



dq 



■\j{Ek{q) - fi n ) - f(E 1 (q) + Mn )], (7.51) 



where f(E) are Fermi functions, and the factor of 2 is 
from the spin sum. Finally, upon differentiating with re- 
spect to \x n and performing the momentum integral one 
extracts the desired result for the low temperature su- 
perfluid stiffness: 



K S (T) = K° s - c—k B T, (7.52) 
with the dimensionless constant c = (ln2/2n). 



VIII. VORTICES 

A. hc/2e versus hc/e vortices 

Having successively incorporated phase fluctuations 
into the effective low energy description of the d x i_ y i 
state, we now turn to a more interesting task - quantum 
disordering the superconductivity to obtain the nodal liq- 
uid phase, a novel Mott insulator. The superconductivity 
is presumed to be destroyed by strong quantum fluctu- 
ations of the order parameter phase tp driven by vor- 
tex excitations. In two-dimensions vortices are simply 
whorls of current swirling around a core region. But in a 
superconductor the circulation of such vortices is quan- 
tized, since upon encircling the core the phase tp can only 
change by integer multiples of 2tt. Inside the core of a 
vortex the magnitude of the complex order parameter | A| 
vanishes, but is essentially constant outside. In the su- 
perconducting phase, the size of the core is set by the 
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coherence length - roughly 10A in the Cuprate materi- 
als. Such vortices are thus tiny "point-like" objects, with 
a truly microscopic size in the Cuprate materials. 

The "elementary" vortex has a phase winding of ±2ir. 
When a superconductor is placed in an external magnetic 
field, the currents circulating around the core of a vortex 
tend to screen out the magnetic field, except within a 
region of the penetration length, A, from the vortex core. 
(In the cuprate materials A is in the range of a thousand 
angstroms.) In addition to the circulation, the total mag- 
netic flux near a vortex is quantized - in units of the flux 
quantum hc/2e. An "elementary" vortex quantizes pre- 
cisely hc/2e of magnetic flux, and will thus henceforth, be 
referred to as an hc/2e vortex. As we shall argue,0-3 to 
obtain the nodal liquid phase it will be necessary to "lib- 
erate" double-strength hc/e vortices, keeping the hc/2e 
vortices "confined". 

Generally, the position of these "point-like" vortices 
can change with time, and their dynamics requires a 
quantum mechanical description. Thus a collection of 
many vortices can be viewed as a many body system of 
"point-like" particles. Since positive (+1) and negative 
(— 1) circulation vortices can annihilate - and disappear 
(just as for real elementary particles like electrons and 
positrons), they behave as "relativistic" particles. There 
is a conserved vortex "charge" in this process, namely 
the total circulation, and an associated current. Since 
the Cooper pairs are Bosons, one anticipates that the 
"dual particles" - the vortices - are also Bosonic forming 
a relativistic Boson system, and this is indeed the case.13 

However, in the superconducting phase at zero temper- 
ature there are no vortices present - this phase constitutes 
a "vacuum" of vortices. More precisely, due to quantum 
fluctuations vortices are present as short-lived "virtual" 
fluctuations, popping out of the "vacuum" in the form of 
small tightly bound (neutral) pairs. For the low energy 
properties of the superconductor these fluctuations can 
be largely ignored. But what happens if these virtual 
pairs unbind into a proliferation of free mobile vortices? 
Vortex motion is very effective at scrambling the phase ip 
of the superconducting order, so that mobile vortices will 
in fact destroy the superconductivity. Since the vortices 
are Bosonic, once they are free and mobile they will "Bose 
condense", at least at zero temperature. One thereby 
obtains a non-superconducting insulating state, with the 
"vortex-condensate" serving as an appropriate order pa- 
rameter. As we shall see, it will be-extremely convenient 
to pass to a "dual" representationE£lc2l in which the vor- 
tices are the basic "particles" - rather than the Cooper 
pairs. 

Consider first unbinding and condensing the "elemen- 
tary" hc/2e vortices.E£l When a Cooper pair is taken 
around such a vortex it's wave function acquires a ±2ir 
phase change. Likewise, when an hc/2e vortex is taken 
around a Cooper pair, the vortex wavefunction acquires 
the 2ir phase change. Thus, hc/2e vortices "see" Cooper 
pairs as a source of "dual flux", each carrying one unit. 
(This notion can be made precise by performing a dual- 



ity transformation - see below and the Appendix.) For 
a Hubbard model of electrons at half-filling, on average 
there is one-half of a Cooper pair per site, as seen ex- 
plicit ly in the effective lattice Cooper pair Hamiltonian, 
Eqn. 7.9, which has offset charge no = 1/2. Thus, these 
elementary vortices "see" a dual "magnetic field", with 
one-half of a dual flux-quantum per plaquette. When the 
hc/2e vortices unbind and condense, they will quantize 
this dual flux, in precisely the same way that the con- 
densation of Cooper pairs in a real superconductor will 
quantize an applied magnetic field - forming an Abrikosov 
flux-lattice (if Type II). The analog of the Abrikosov 
flux-lattice for the hc/2e vortex condensate is an ordered 
lattice of Cooper pairs. In this "crystal" state at half- 
filling, the Cooper pairs will preferentially sit on one of 
the two equivalent sub-lattices of the square lattice. This 
state can be described as a commensurate charge-density- 
wave with ordering wavevector (Q = tt, 7r), which sponta- 
neously breaks the discrete symmetry under translation 
by one lattice spacing. Such ordering implies a consid- 
erable degree of double occupancy for the electrons, and 
thus seems most reasonable for a Hubbard type model 
with an attractive on-site interaction (negative u) . In the 
Cuprate materials there is a strong on-site repulsion, and 
moreover there is no evidence for "charge-ordering" near 
Q. Thus, for a description of the pseudo-gap regime in 
the Cuprate materials, we can rule out the hc/2e vortex- 
condensate on phenomenological grounds. 

Instead, we consider the possibility of unbinding and 
condensing double-strength hc/e vortices, keeping the el- 
ementary hc/2e vortices confinedJlj When an hc/e vor- 
tex is taken around a Cooper pair it acquires a 4-7T phase 
change. A 2ir phase change corresponds to taking such 
an hc/e vortex around "half of a Cooper pair - which has 
charge e. Thus, a condensation of hc/e vortices should 
correspond to a "crystal" of such charge e objects. But 
at half-filling with charge e per lattice site, this should 
correspond to a state without charge ordering or transla- 
tional symmetry breaking. As we shall see, for a d x 2_ y 2 
superconductor the resulting hc/e "vortex-condensate" 
gives a description of the nodal liquid phase. 

This procedure - keeping the elementary hc/2e vor- 
tices confined and only liberating the hc/e vortices - is 
responsible for the remarkable properties of the nodal 
liquidJl3Tb see why, consider first the Berry's phase term 
in Eqn. 7.3. With only hc/e vortices present, the Cooper 
pair phase, ip, only winds by integer multiples of 47r - 
not 2tt. At half-filling (with no = 1/2) the Berry's phase 
term will not contribute to the partition function (see 
Eqn. 7.7) and can thus be dropped entirely in the de- 
scription of the nodal liquid. This can be implemented 
by defining a new phase field: 



= p/2, 



(8.1) 



and only allowing vortices in 4>{x) with circulation 2n 
times an integer. This restriction precludes hc/2e vor- 
tices, and guarantees that the field 
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J<t> 



is single-valued. As an operator, b creates a spinless exci- 
tation with charge e. When re- written in terms of </>, the 
effective Lagrangian for a d-wave superc onductor with 
quantum phase fluctuations (from Eqns. 7.35 , 7.36] ) be- 
comes L = La, + Li n t + L-j, with 



Li 



-K u (d a <j)- A a 



(8.3) 



and Lip given in Eq. 7.28| . The Berry's phase term has 
been dropped, since it plays no role when exp{i<j)) is a sin- 
gle valued field. Here, we have absorbed a factor of two 
into n u and also completed the square with the nodon 
current, J a , dropping order terms which are irrelevant 
as discussed after Eq. 7.43 . Notice that the coefficient of 
A a is one - as expected for a charge e operator exp(iip). 
By precluding hc/2e vortices, we see the emergence of 
a new bosonic field, exp(i(j)), with exotic quantum num- 
bers - charge e but spin zero - which will be referred 
to as a "holea". This is the first hint of spin-charge 
separationeirEil in the nodal liquid. 

As we shall see, another remarkable consequence of 
precluding hc/2e vortices, is that the charge neutral spin 
one-half nodons survive under hc/e vortex condensation 
into the nodal liquid. To see why this is jupt the case 
if elementary hc/2e vortices are condensed!!] (as in the 
charge-density- wave) , it is very instructive to consider 
the transformation which relates the nodons to the d- 
wave quasiparticles, Eq. |7.26] , which can be written in 
terms of the new field 4> (= <p/2) as: 



(8.4) 



In the presence of vortices, the nodon field ip only re- 
mains single- valued if hc/2e vortices are excluded (so 
that exp(±i(j)) is single valued). Indeed, when a nodon is 
taken around an hc/2e vortex, it's wavefunction changes 
sign, since (f> winds by tt. This implies a very strong and 
long-ranged "statistical" interaction between nodons and 
hc/2e vortices. If hc/2e vortices proliferate and condense, 
it will clearly be very difficult for the nodons-tp propogate 
coherently. In fact, we have argued recentlytij that in this 
case the nodons are bound (actually "confined") to the 
holons, leaving only the electron in the spectrum of the 
charge-density- wave. 



B. Duality 



We now consider implementing the procedure of un- 
binding and condensing hc/e vortices in the d x 2_ y 2 su- 



perconductor. To this end, it is extremely convenient to 
pass to the "dual" representationEjo in which the vor- 
tices are the basic "particles", rather than the Cooper 
pairs. The most straightforward way to incorporate hc/e 
vortices is by placina-the (single-valued) field exp(i(f>) on 
the sites of a latticeB3 so that vortices can exist in the 



plaquettes. A lattice duality transformation can be imple- 
mented in which the phase <p is replaced by a dual field, 
9, which is the phase of a vortex complex field, $ ~ e lB '. 
In a Hamiltonian description, <E> and & can be viewed 
as vortex quantum field operators - which destroy and 
create hc/e vortices. On a 2 + 1-dimcnsional Euclidian 
space-time lattice, the appropriate model corresponding 
to the phase Lagrangian Eq. 3.3 is essentially a classical 
3d-xy model with an effective gauge field: 



A e ff — A — 



(8.5) 



The lattice duality transformation for the 3d-xy model 
with gauge field is implemented in some detail in the 
Appendix. An alternative method which we sketch be- 
low, involves implementing the duality transformation di- 
rectly in the continuum. Cj 

To this end we introduce a vortex 3-current, j^, which 
satisfies, 



J i 



d v d\ 



(8.6) 



In the presence of hc/e vortices, <p is multi- valued, d u <j) is 
not curl-free, and f a is non-vanishing. Even in the dual 
vortex representation the total electrical charge must be 
conserved. This can be achieved by expressing the total 
electrical 3-current (in units of the electron charge e) as 
a curl, 



jtot Q 



•7) 



where we have introduced a "fictitious" dynamical gauge 
field, a u . (In the Appendix the electrical 3-current is ex- 



pressed as a l atti ce curl of a^.) Upon combining Eqn. 7.38 
with |S.6| and ^Tt], one can eliminate the phase field, 4>, and 
relate a a to the vortices: 

fa = tfivxdviK^EXapdaap + A\ K~ X J A ], (8.8) 

where J u is the quasiparticle 3-current defined earlier in 



Eqs.|L3§|73l 



In this continuum approach to duality, a dual descrip- 
tion is obtained by constructing a Lagrangian, Ld, de- 
pending on dju, J a and j™, whose equation of motion, 
obtained by differentiating the action with respect to a M , 
leads to the above equation. It is convenient to first ex- 
press the vortex 3-current in terms of a complex field, 
which can be viewed as an hc/e vortex destruction op- 
erator. The dual Lagrangian is constructed to have an 
an associated U(l) invariance under $ — > e la Q, which 
guarantees that f is indeed conserved. When an hc/e 
vortex is taken around a Cooper pair it aquires a 47r phase 
change (2tt around a charge e "holon" ) . In the dual repre- 
sentation the vortex wavefunction $ should acquire a 4tt 
phase change (or 2tt for a "holon"). This can be achieved 
by minimally coupling derivatives af $ to the "fictitious" 
vector potential a u . 

The appropriate dual Lagrangian can be conveniently 
decomposed as Ld = L^, + L v + L a , where L^ is given 
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in Eq 47.2S . The vortex piece has the Ginzburg-Landau 
form,0 



C v = !f\(d li -iaJ*\ ,, -V*m), 



(8.9) 



as constructed explicitly with lattice duality in the Ap- 
pendix. The vortex 3-current, following from — 
-dCv/dcLft, is 



3l = K M Im[**(a^ 



(8.10) 



For small |$| (appropriate close to a second order transi- 
tion) one can expand the potential as, V$(X) = r$X 2 + 
u$X . The remaining piece of the dual Lagrangian is 

C a = — (e 2 j - b 2 ) + a^xduiAx - a^ 1 J A ), (8.11) 

with dual "magnetic" and "electric" fields: b = eijdiaj 
and ej — v^r 1 (djao — d$aj). It can be verified that the 



dual Lagrangian has the desired property that Eq. 8.8 



follows from the equation of motion SSo/da^ = 0. 



IX. NODAL LIQUID PHASE 



In this section we employ the dual representation of 



the d 



superconductor to analyze the quantum dis- 
ordered phase - the nodal liquid. The dual representation 
comprises a complex vortex field, which is minimally cou- 
pled to a gauge field, as well as a set of neutral nodon 
fermions. Without the nodons and in imaginary time, 
the dual Lagrangian is formally equivalent to a classical 
three-dimensional superconductor at finite temperature, 
coupled to a fluctuating electromagnetic field. To disor- 
der the d-wave superconductor, we must order the dual 
"superconductor" - that is, condense the hc/e vortices. 
The nature of the resulting phase will depend sensitively 
on doping, since upon doping, the dual "superconductor" 
starts seeing an applied "magnetic field" . Below, we first 
consider the simpler case of half-filling. We then turn to 
the doped case, where two scenarios are possible depend- 
ing on whether the dual "superconductor" is Type I or 
Type 1IM 



A. Half-filling 

Specialize first to the case of electrons at half-filling, 
with particle-hole symmetry. In the dual representation, 
the "magnetic field" , b, is equal to the deviation of the to- 
tal electron density from half-filling. Thus at half-filling 
(b) = and the dual Ginzburg-Landau theory is in zero 
applied field. The quantum disordered phase corresponds 
to condensing the hc/e vortices, setting ($) = $ 7^ 0. In 
this dual Meissner phase the vortex Lagrangian becomes 



(9.1) 



where a* represents the transverse piece of a^. It is then 
possible to integrate out the field a M which now enters 
quadratically in the Lagrangian. Equvalently, a M can be 
eliminated using the equation of motion which follows 
from SSd /Sa^ = 0. The full Lagrangian in the nodal 
liquid phase is then 



2^o 



0[{dJf 



(9.2) 



where we have introduced the physical magnetic and elec- 
tric fields: B — eijdiAj and Ej = djAg — dtAj. The last 
two terms describe a dielectric, with magnetic perme- 
ability /io = Ko'J'o an d dielectric constant eo = (/We) -1 ; 
with the sound velocity entering, rather than the speed 
of light. The external electromagnetic field is coupled to 
the 3-current 1^ , which can be expressed as a bi-linear of 
the nodon fermions as, 



(9.3) 



^0 v c 



Notice that this 3-current is automatically conserved: 
V M = 0. 

The order (dJ) 2 terms which we have not written out 
explicitly are quartic in the fermion fields, and also in- 
volve two derivatives. Since describes Dirac fermions 
in 2 + 1 space-time dimensions, these quartic fermion 
terms arc highly irrelevant, and rapidly vanish under the 
rescaling transformation in Eqn. 7.43. Thus, the low en- 



ergy description of the nodal liquid phase is exceedingly 
simple. It consists of four neutral Dirac fermion fields - 
two spin polarizations (a = 1, 2) for each of the two pairs 
of nodes. Despite the free fermion description, the nodal 
liquid phase is highly non-trivial when re-expressed in 
terms of the underlying electron operators. Indeed, the 
ij) fermion operators are built from the quasiparticle oper- 
ators ^ in the d-wave superconductor, but are electrically 
neutral, due to the "gauge transformation" in Eq. 7.26. 



In addition to the gapless nodons, one expects exotic 
charged excitations at finite energy in the nodal liquid. 
To see this, imagine applying an external dual "mag- 
netic field" to the Ginzburg-Landau "superconductor", 
which corresponds to a non-zero chemical potential for 
the electrons. Being in the Meissner state, this "field" 
will be screened out, so that the internal field, b, which 
corresponds to deviations in the electron charge density 
from half-filling, j-«fill vanish. Clearly, this corresponds to 
a Mott insulatonla with the Mott gap being proportional 
to the dual critical field. In a Type II superconductor, an 
internal magnetic field will beJ'quantized" into flux-tubes 
carrying a quantum of fluxo. For the dual Ginzburg- 
Landau theory, this corresponds to a quantization of elec- 
tric charge, with a flux tube corresponding to charge e. 
Thus, in the nodal liquid one expects the presence of 
gapped finite energy excitations with charge e. These 
"holon" excitations are exotic since they carry no spin. 
The holon is the basic topological excitation that can be 
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created in the hc/e vortex-condensate. The existence of 
a spin one-half neutral nodon excitation and a spinless 
charge e holon excitation in the nodal liquid, _is .-a, dra- 
matic demonstration of spin-charge separation.circfj The 
excitations in the nodal liquid have the same quantum 
numbers as in the spin-charge separated gauge theories ,e!j 
but are weakly interacting, rather than strongly coupled 
by a gauge field. 



1. Spin response 

Although the nodons are electrically neutral they do 
carry spin, so the low-energy spin response in the nodal 
liquid can be computed from the Dirac Lagrangian C^. 
Moreover, since C$ was not altered under the duality 
transformation, the spin properties of the nodal liquid 
are essentially identical to those in the d x i_ y i supercon- 
ducting phase. As a simple example, consider the uni- 
form magnetic spin susceptibility, x- The uniform part 
of the electron spin oper ator is given as the conserved 



spin density in Eqn. 7.41 



1 



(9.4) 



Being bi-linear in nodon operators spin correlation and 
response functions can be readily computed from the free 
nodon theory For example, the uniform spin suscepti- 
bility is given by 



X 



dE(-df/dE)p n (E), 



(9.5) 



where the nodon density of states is p n (E) = 
(const) E/vfV a , and f(E) is a Fermi function. One finds 
X ~ T/vpv^. There are also low energy spin excita- 
tions at wavevectors which span between two different 
nodes. The associated spin operators can be obtained by 
re-expressing the electron spin operator, 



-fc! 



(9.6) 



in terms of the nodons. For example, the staggered mag- 
netization operator, S^, is found to be 



s ^ = \ [$(T v tro v )ift +h.c] 



(9.7) 



Notice that this operator is actually "anomalous" in 
terms of the conserved nodon charge. 

In addition to carrying spin, the nodons carry energy, 
and so will contribute to the thermal transport. In the 
absence of scattering processes (such as Umklapp) the fi- 
nite temperature nodon thermal conductivity is infinite. 
In practive, impurities will scatter the nodons and lead 
to a finite thermal conductivity. In fact, impurity scat- 
tering should also play an important role in modifying 
the spin response of the nodal liquid. 



2. Charge response 



The electrical charge properties of the nodal liquid are 
of course very different than in the superconductor. To 
see this, imagine changing the chemical potential away 
from p = which corresponds to applying an exter- 
nal "magnetic" field to the dual Ginzburg-Landa theory: 
£ M = —fib. Being in the "Meissncr" phase, the elec- 
tron density will stay "pinned" at half-filling for p < p c , 
with fi c the Ginzburg-Landau critical field. Despite the 
presence of this charge gap, there are low energy current 
fluctuations in the nodal liquid. Indeed, in this phase 
the electrical current operator is 1^, which is bi-linear 
in the nodon fermions, if). To compute the electrical 
conductivity in the nodal liquid requires computing a 
two-point correlator of J M at zero wavevector (say in the 
x— direction) I x (q = 0) = (e /K o v 2 )<9 2 J x (q — 0). But no- 
tice that J x (q = 0) is proportional to a globally conserved 
nodon charge, since J x (x) = (vp /2)t(j 1 ij} 1 . Thus, when 
the nodon number is conserved one has I x (q = 0) = 0, 
and the nodons do not contribute to the electrical con- 
ductivity. When impurity (or Umklapp) scattering is 
present, however, the nodon number is no longer con- 
served, and the nodons will contribute to the real part of 
the electrical conductivity, but only at finite frequencies. 

It is instructive to briefly consider the behavior of 
the electron Green's function, which can be accessed in 
photo-emission and tunneling experiments. The electron 
operator c a (x) can be decomposed as a product of nodon 
and holon operators. For example, near the node at K j 
one can write, 



c a (x) 



(9.8) 



where ip is a nodon operator and exp(i<f>) can be inter- 
preted as a holon destruction operator. In the nodal 
liquid phase, the electron Green's function, G(x,t) = 
(c\x,t)c(0,0)) factorizes as, 



G(x,t) 



iKi 



i^,t) e mm)G n {x,t), (9.9) 



where the nodon Green's function is, 



(9.10) 



Although G n (x,t) decays as a power law |ir|~ 2 and i -2 , 
since creating a holon costs a finite energy the holon 
Green's function is expected to be short-ranged, decay- 
ing exponentially in space and time. This indicates a gap 
in the electron spectral function at the Fermi energy. 
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B. Doping the Nodal Liquid 

We briefly discuss the effects of doping charge into the 
nodal liquid phase. In a grand canonical ensemble this is 
achieved by changing the chemical potential, fi = Aq. In 
the dual Ginzburg-Landau description of the vortices, a 
chemical potential acts as an applied dual field, as seen 
from Eq. |S.llj 



since 



-[lb. 



(9.11) 



The dual magnetic field, b = e^diCLj, is the total elec- 
tric charge in units of e. Provided the applied dual 
field, fi, is smaller than the critical field (/x c ) of the 
Ginzburg-Landau theory, the dual superconductor stays 
in the Meissner phase - which is the nodal liquid phase 
at half-filling. But for \i > fj, c dual flux will penetrate the 
Ginzburg-Landau superconductor, which corresponds to 
doping the nodal liquid. The form of the dual flux 
penetration will depend critically on whether the dual 
Ginzburg-Landau theory is Type I or Type II. Within a 
mean-field treatment this is determined by the ratio of 
the dual penetration length, A„, to the dual coherence 
length, (where the subscript v denotes vortices). In 
particular, Type II behavior is expected if \ v /£ v > l/v2, 
and Type I behavior otherwise. In the Ginzburg-Landau 
description A„ determines the size of a dual flux tube, 
which is essentially the size of a Cooper pair. We thus 
expect that X v will be roughly equal to the supercon- 
ducting coherence length, £, which is perhaps 10 — 15A 
in the cuprates. On the other hand, is the size of the 
"vortex-core" in the dual vortex field, and presumably 
can be no smaller than the microscopic crystal lattice 
spacing, £ v > 3 — 5 A. This reasoning suggest that \ v /£ v 
is probably close to unity in the cuprates, so that either 
Type I or Type II behavior might be possible - and could 
be material dependent. We first consider such Type II 
doping, returning below to the case of a Type I Ginzburg- 
Landau theory. 



1. Type II Behavior 

The phase diagram of a clean threes/dimensional type 
II superconductor is well understood.O Above the lower 
critical field, H e ±, flux tubes penetrate, and form an 
Abrikosov flux lattice - usually triangular. As the applied 
field increases the flux tubes start overlapping, when their 
separation is closer than the penetration length. Upon 
approaching the upper critical field H C 2 their cores start 
overlapping, the Abrikosov flux lattice disappears, and 
the superconductivity is destroyed. 

These results hold equally well for our dual Ginzburg- 
Landau superconductor, except that now the direction 
parallel to the applied field is actually imaginary time. 
Moreover, the Ginzburg-Landau order parameter de- 
scribes quantum (hc/e) vortices, and the penetrating flux 
tubes are spinless charge e holons. Upon doping the 



nodal liquid with /i > fj, c ±, charge is added to the 2d 
system, which corresponds to the penetration of dual 
magnetic flux. In this dual transcription, the resulting 
Abrikosov flux-lattice phase is a Wigner crystal of holons, 
with one holon per real space unit cell of the lattice. 
Upon further doping, at fi — /i C 2, the crystal of holons 
melts, and they condense - this is the d-wave supercon- 
ductor. 

In the holon Wigner crystal phase, translational sym- 
metry is spontaneously broken. However, in a real mate- 
rial the Wigner crystal will have a preferred location, de- 
termined by impurities and perhaps crystal fields, which 
will tend to pin and immobilize the crystal. The resulting 
phase should be an electrical insulator. 

A striking and unusual feature of the holon Wigner 
crystal is that it co-exists with the nodal liquid. We 
thereby arrive at a description of a rather remarkable 
new phase of matter. A Wigner crystal of doped holons 
co-exists with neutral gapless fermionic excitations - the 
nodons. In this co-existing phase, low energy spin and 
thermal properties will be dominated by the nodons. The 
behavior will be qualitatively similar to that in the un- 
doped nodal liquid phase. It is possible that this phase 
underlies the physics of the pseudo-gap region of the high 
T c cuprates. 



2. Type I Behavior 

In a Type I superconductor, the applied field is ex- 
pelled until the critical H c is exceeded. a At this point 
there is a first order phase transition from the Meissner 
phase with all the flux expelled, to a normal metal phase 
in which (essentially) all the field penetrates. If our dual 
Ginzburg-Landau theory is type I, then analogous prop- 
erties are expected. Specifically, as the chemical poten- 
tial increases, the dual field - which is the holon density 
- remains at zero until a critical chemical potential fi c is 
reached. At this point there is a first order phase transi- 
tion, between the nodal liquid phase at half-filling, and a 
d-wave superconductor at finite doping, x c . At fixed dop- 
ing x < x c , phase separation is impeded by long-ranged 
Coulomb interactions between the holons. The system 
will break apart into co-existing "micro-phases" of nodal 
liquid and d-wave superconductor. The configuration of 
the "micro-phases" will be determined by a complicated 
competition between the Coulomb energy and the (posi- 
tive) energy of the domain walls. In practice, impurities 
will also probably play a very important role. 



C. Closing Remarks 

The theoretical framework described above gives a 
skeletal description of the nodal liquid and, upon doping, 
the holon Wigner crystal. There are many important is- 
sues which will need to be addressed in detail to sec if this 



27 



novel Mott insulating phase gives a correct description of 
the low temperature pseudo-gap regime in the cuprates. 
At very low doping the cuprates are antiferromagnetic so 
it will clearly be necessary to incorporate magnetism into 
the theoretical framework. Perhaps even more important 
is assessing the role of impurities, which are expected to 
have rather dramatic effects both on the holon Wigner 
crystal and the gapless nodons. Impurities will tend to 
disorder the Wigner crystal and will scatter the nodons 
probably leading to a finite density of states and diffusive 
rather than ballistic motion. Since the nodons carry spin 
but no charge, a rather exotic "spin metal" phase is pos- 
sible with a finite "spin conductivity" (but zero electrical 
conductivity) even at zero temperature. It is also possible 
that the impurities will localize the nodons, perhaps lead- 
ing to a random singlet phase or a spin glass. An addi- 
tional complication is that some materials might exhibit 
phase separation upon doping (Type I rather than Type 
II behavior) exhibiting micro-phase co-existence between 
the antiferromagnet and the d-wave superconductor, pre- 
empting the nodal liquid phase. It clearly remains as a 
future challenge to fully sort out the mysteries of the 
pseudo-gap regime. 

A more general theme of these notes is that novel spin 
liquid phases can sometimes be more conveniently viewed 
as descendents of superconductors - rather than the more 
traditional route via magnetism. One can imagine quan- 
tum disordering other exotic superconducting phases be- 



sides the dr. 



state, to obtain new spin liquid phases. 



Perhaps some of these phases will appear in other systems 
which exhibit finite angular momentum pairing, such as 
3 — He and the heavy Fermion materials. 
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APPENDIX A: LATTICE DUALITY 

Duality plays a key role in understanding how to quan- 
tum disorder a superconductor, both in 1 + 1 space-time 
dimensions (Sec. 0) and in 2 + 1 (Sec. |VIlH ). The key 
idea involves exchanging the order parameter phase </> for 
vortex degrees of freedom. In r L + 1 dimensions these are 
point- like space-time vortices, La whereas in 2 + 1 then&^ire 
poin t like vortices in space which propogate in time.EJ In 
Sec. VIII we chose to work directly in the continuum in 
implementing the 2 + 1 duality transformation. However, 
the physics of duality is perhaps more accessible when 
carried out on the lattice. In this Appendix we show in 
some detail how lattice duality is implemented in both 
1 + 1 and 2 + 1 dimensions .tilo For simplicity we first 
Wick rotate to Euclidian space, and rescale imaginary 
time to set the charge velocity to one. The appropriate 
lattice model is then simply a 2d square lattice or 3d cu- 
bic lattice xy model. In the latter case, we also want to 
include a gauge-field, A, which is a sum of the physical 
electromagnetic field and the nodon current, as discussed 
in Sec. - see Eq. pi 



The degrees of freedom which live on the sites of the 
square or cubic lattice (denoted by a vector of integers 
x) are the phases 4> x G [0, 2w]. As usual, the gauge field 
lives on the links. Discrete lattice derivatives are denoted 
by 

A^x = (f> x+fl - (f> x , (Al) 

where \x = x,y for the square lattice and fj, — x,y, z for 
the cubic lattice and x + \i denotes the nearest neighbor 
site to x in the fi, direction. The gauge field is minimally 
coupled via, 

A^-* A^ + A*. (A2) 



Consider the partition function, 



Z= / TT d<f) x exp[y2 V K (A M ^ X )]. 
Jo 



(A3) 



Here the periodic "Villain" potential V K is given by, 

oo 

esp^A^)] = J2 e- Kj2 / 2 e UA *, (A4) 



J— — oo 

with integer J. When k >> 1 only the terms with 
J = 0, ±1 contribute appreciably in the sum and this 
reduces to the more familiar form: 

V K (A4>) = Kcos(Act>), (A5) 

with K = 2exp(-n/2). 

The partition function can thus be expressed as a sum 
over both (f> and a vector of integers, J Xl with components 
J£ living on the links of the lattice: 



z= [ n^E e_5 ^ Tr 



4>,J e 



(A6) 



with action 
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S = S a + ^2i(K-f ai )<l> x , (A7) 

X 

^o-fEi^i 2 - ( A8 ) 

x 

In this form the integration over <j> can be explicitly per- 
formed giving 

Z = Tr' f e- So , (A9) 

where the prime on the trace indicates a divergenceless 
constraint at each site of the lattice: 

A • J x = 0. (A10) 

In the presence of a gauge field there is an additional 
term in the action of the form, 

S A = iY,Jx-A x . (All) 

X 

It is thus clear that the integer of vectors J can be in- 
terpreted as a conserved electrical current flowing on the 
links of the lattice. The divergenceless constraint on this 
electrical 3-current can be imposed automatically by re- 
expressing J as a curl of an appropriate dual field. Con- 
sider first the 2d case. 



1. Two dimensions 

To guarantee divergenceless we set the current equal 
to the (2d) curl of a scalar field, 9 X : 

2ttJ£ = e^A„9 x , (A12) 

so that the action becomes 

S°W = ^E( A A) 2 - (A13) 

X,fJ, 

To insure that J is an integer field, 9 must be constrained 
to be 27r times an integer. This additional constraint can 
be imposed by introduction of yet another integer field, 
n x , which will be interpreted as the (space-time) vortex 
density. The partition is thereby re-expressed as (drop- 
ping an unimportant multiplicative constant), 




(A14) 



with 

S = S {e) + Y\~nl+in x 6 x \. (A15) 

X 

For k = the summation over n x gives a sum of delta 
functions restricting 9 x /2tt to be integer. But we have 



softened this constraint, introducing a vortex "core" en- 
ergy K ^ 0. 

At this stage one could perform the Gaussian integral 
over 9, to obtain a logarithmically interacting plasma of 
(space-time) vortices. Alternatively, for h » 1 the sum- 
mation over n x can be performed giving, 

S = S Q (9)-uY,cos(9 x ), (A16) 

X 

with u = 2exp(—n/2). Upon taking the continuum limit, 
9 x — > 9(x), one recovers the (Euclidian) sine-Gordon the- 
ory, S — J d 2 xC with 

£ = -^(V9) 2 -ucos(9). (A17) 

After Wick rotating back to real time and restoring the 
velocity this takes the identical form to the dual La- 
grangian considered for the 2-leg ladder in Sec. |v|. 

2. Three dimensions 

In three dimensions the divergenceless integer 3- 
current J can be written as the curl of a vector field, 
a: 

2ttJ x = A x a x . (A18) 

As in 2d one imposes the integer constraint (softly) by 
introducing an integer vortex field, in this case a 3- vector 
j, to express the partition function as, 

/oo 
J|da^e- S , (A19) 

with 

S = S a (a) +Y / [^\3x\ 2 ~ ijx ■ S x ], (A20) 



So(a) = £^Y,\ Kx ^\ 2 - (A21) 

X 

The integer vector field j is the vortex 3-current, "mini- 
mally" coupled to a. To see that the vortex 3-current is 
conserved, it is convenient to decompose the vector field 
a into transverse and longitudinal pieces: a = at — A#, 
with 9 X a scalar field. The action becomes, 

S = S (a) + ^\j x \ 2 + ij x ■ (A9 X ~ a x )\, (A22) 

X 

where we have dropped the subscript "t" on a. The par- 
tition function follows from integrating over both a and 
9 and summing over integer j. Integrating over 9 leads 
to the expected condition: A • j = 0. Alternatively, for 
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k » 1 one can perform the summation over j to arrive 
at an action depending on 6 and a: 

5 = 5 (a)-A'^cos(A,A-<), (A23) 

X,fJ- 

with K = 2exp(-R/2). 

In the presence of a gauge field there is an addi- 
tional term in the action of the form, 



S A 



i 

2^ 



J^(A x a x ) ■ A x , 



(A24) 



which follows directly from Eqn. [All| and [AlS . 

At this stage one can take the continuum limit, let- 
ting a x — > a(x) and X — > 9(x). Upon expanding the 
cosine for small argument one obtains S = J d 3 x£ with 
(Euclidian) Lagrangian 



£ = ^(Vxa) 2 + ^(W-a) 



(A25) 



In this dual representation, the vortex 3-current (which 
follows from dC/da) is given by j v = K(W8 — a). Notice 
that the vortices are minimally coupled to the "vector 
potential" a, whose curl equals the electrical 3-current. 
The field 9 can be interpreted as the phase of a vortex 
operator. In fact it is convenient to introduce such a 
complex vortex field before taking the continuum limit: 



(A26) 



The continuum limit can then be taking retaining the full 
periodicity of the cosine potential. The appropriate vor- 
tex Lagrangian replacing the second term in Eqn. A25 
is, 



£, = y|(V-^)$| 2 + Vi(|$|). 
The vortex current operator becomes, 



J 



A7m["T(V -ia)®]. 



(A27) 



(A28) 



If the potential is expanded for small $ as V$(X) — 
r$A 2 + m$A 4 , the full dual theory is equivalent to a 
Ginzburg-Landau theory for a classical three-dimensional 
superconductor. Inclusion of the original gauge field 
leads to an additional term in the dual Lagrangian: 



£a = t-(V x a) ■ A. 
Zn 



(A29) 



After Wick rotating back to real time and restoring the 
velocity, C + Ca be com es i denti cal to the dual vortex 



Lagrangian in Eqns. 8.S and 3.11 
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